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O \ Abstract 

■ Heavy neutral Higgs boson production and decay into neutralino and chargino pairs 

is studied at the Large Hadron Collider in the context of the minimal supersymmetric 
standard model. Higgs boson decays into the heavier neutralino and chargino states, 
W\ ■ i-^-^ A^ ^ x'ixl Xlxl Xlxl Xlxl xlxl as well as , A^ ^ xfx^,xtx2 (all lead- 



ing to four-lepton plus missing transverse energy final states), is found to improve the 
possibilities of discovering such Higgs states beyond those previously identified by con- 
sidering H^,AP X2X2 decays only. In particular, H^,AP bosons with quite heavy 
masses, approaching ~800 GeV in the so-called 'decoupling region' where no clear SM 
signatures for the heavier MSSM Higgs bosons are known to exist, can now be discerned, 
for suitable but not particularly restrictive configurations of the low energy supersym- 
metric parameters. The high Ma discovery reach for the and A^ may thus be greatly 
extended. Full event-generator level simulations, including realistic detector effects and 
analyses of all significant backgrounds, are performed to delineate the potential H^,J^ 
discovery regions. The wedgebox plot technique is also utilized to further analyze the 
plus missing transverse energy signal and background events. This study marks the first 
thorough and reasonably complete analysis of this important class of MSSM Higgs boson 
signature modes. In fact, this is the first time discovery regions including all possible 
neutralino and chargino decay modes of the Higgs bosons have ever been mapped out. 



1 Introduction 



Among the most investigated extensions of the standard model (SM) are those incorporating 
supersymmetry (SUSY), and among these the one with the fewest allowable number of new 
particles and interactions, the minimal supersymmetric standard model (MSSM), has certainly 
received considerable attention. Yet, when prospective signals at the Large Hadron CoUider 
(LHC) of the new particle states within the MSSM are considered, there is still much that 
needs clarification. Nothing underscores this more than the MSSM electroweak symmetry 
breaking (EWSB) Higgs sector. Included therein is a quintet of Higgs bosons left from the 
two SU{2)^ Higgs doublets after EWSB (see [1, 2] for more details): a charged pair, H^, 
the neutral CP-odd and the neutral CP-even and (with Mh < Mh)- The entire 
Higgs sector (i.e., masses and couplings to ordinary matter) can be described at tree-level by 
only two independent parameters: the mass of one of the five Higgs states {e.g.^ Ma) and the 
ratio of the vacuum expectation values of the two Higgs doublets (denoted by tan/5). These 
must be augmented to include significant radiative corrections which most notably raise the 
upper limit on the mass of the light Higgs boson from < Mz at tree-level to < 140 GeV 
(150 GeV) with inclusion of corrections up to two loops and assuming a stop-sector scale of 
MsusY = ITeV (2TcV) and rrit = (178.0 ± 4.3) GeV according to [3], or <135 GeV with 
rrit = (172.6 ± 1.4) GeV by [4] (stop mass range not specified). This definite upper bound will 
allow experimentalists to definitively rule out such a minimal SUSY scenario at the LHC if 
such a light Higgs state is not observed. Thus, the possible production and decay modes of 
the /i° state have understandably been investigated in quite some detail [2]. In contrast, the 
possibilities for the other heavier neutral MSSM Higgs bosons have not been so thoroughly 
examined. Yet it is crucial that the avenues for discovery of these other MSSM Higgs bosons 
be well understood since, even if a candidate for discovery is experimentally identified, it 
may be indistinguishable from a SM Higgs boson (this corresponds to the so-called 'decoupling 
region', with Mh,Ma ^ 200 GeV and for intermediate to large values of tan/3 [2, 5]). Then 
the additional identification of heavier Higgs bosons may well be required to establish that 
there is in fact an extended Higgs sector beyond the single doublet predicted by the SM. 

Finding signatures for these heavier MSSM Higgs bosons has proved to be challenging. 
Unlike the lone Higgs boson of the SM of similar mass, couplings of these MSSM Higgs bosons 
to SM gauge bosons are either absent at tree level (for A^) or strongly suppressed over much of 
the allowed parameter space (for H^). Thus, identification of A^ and via their decays into 
known SM particles relies chiefiy on decays of said Higgs bosons into the heaviest fermions 
available, namely, tau leptons and bottom quarks^. Identification of hadronic decays/jet 
showers of these third generation fermions may be problematic in the QCD-rich environment 
of the LHC^, so that it is very questionable that the entire parameter space can be covered 
with just SM-like signatures. Fortunately, in the MSSM there is an alternative: decays of these 
Higgs bosons into sparticles, in particular the charginos and neutralinos^ formed from the EW 
gauginos and Higgsinos. Higgs boson couplings to certain -ino states may be substantial, and 
these heavy sparticles may themselves decay — except for x\ which is assumed to be the stable 

AP top quark couplings are suppressed relative to a SM Higgs boson of the same mass. 
^In addition, jet-free events from Higgs boson decays to tau-lepton pairs where both tau-leptons in turn 
decay leptonically also come with considerable background-separation challenges [6] . 

^In the remainder, charginos and neutralinos collectively will be abbreviated by '-inos'. 
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lightest supersymmetric particle (LSP) — in readily-identifiable ways (such as into leptons) 
to provide a clean experimental signature. 

A number of previous articles [7, 8, 9, 10, 11] as well as at least one Ph.D. thesis [12] have 
focused on the signal potential of the decays of the heavier neutral MSSM Higgs bosons into 
neutralinos and charginos: 

H',A'^x:x;,X-x' {a,b^l,2, i, j = 1, 2, 3, 4). (1) 

Therein only subsequent -ino decays into leptons (which will be taken to mean electrons 
and/or muons, £ = e,n) were considered, as this is preferable from the standpoint of LHC 
detection. Since relatively light sleptons can greatly enhance [13, 20, 21] the branching ratios 
(BRs) for such decays, the properties of the slepton sector of the MSSM also need to be 
specified. All of the previous works concentrated almost^ exclusively on the decays H^, — 
X2X2- addition, the subsequent ncutralino decays X2 ~^ Xi^"*"^" were typically presumed to 
proceed via three-body decays with an ofT-mass-shell intermediate Z^* or slepton, neglecting 
the possibility of the intermediate or slepton being on-mass-sheU ([22] and [23] delve in 
considerable depth into the distinctions between these cases). 

In this work^, all the decays in (1) are incorporated. In fact, as the presumed mass of 
a Higgs boson grows, more such decay modes will become accessible. Therefore, if decay 
channels to the heavier -inos are significant, they may provide signatures for heavier neutral 
Higgs bosons (with masses well into the aforementioned decoupling region). When heavier 
-ino states are included, it also becomes easier to construct model spectra with slepton masses 
lying below those of the heavier -inos. Thus, in this work, intermediate sleptons are allowed 
to be both on- and off-mass-shell (same for the Z^'^*^)^. More background channels are also 
emulated than in previous studies. The Higgs boson production modes considered herein are 
gg if", (gluon-fusion) and qq — > A^ (quark- fusion) . (The second mode is dominated 
by the case q — h.) 

This work is organized as follows. The next section provides an overview of the MSSM 
parameter space through calculation of inclusive rates for the relevant production and de- 
cay processes contributing to the signal. Sect. 3 then specializes these results to the more 
restrictive minimal supergravity (mSUGRA) scenario for SUSY breaking. Sect. 4 gives the 
numerical results for the signal and background processes based upon Monte Carlo (MC) sim- 
ulations of parton shower (PS) and hadronization as well as detector effects. This includes 
mapping out discovery regions for the LHC. The recently-introduced 'wedgebox' method of 
[27], which is reminiscent of the time- honored Dalitz plot technique, is utilized in Sect. 5 to 
extract information about the -ino mass spectra and the -ino couplings to the Higgs bosons. 
Finally, the last section presents conclusions which can be drawn from this study. 

^The decays H^,AP xtxi 1X1X2 were also studied in [7] but found to be unproductive due to large 
backgrounds to the resulting di-lepton signals. 

preliminary account of this analysis is given in Ref. [24]. 

^Similar studies for charged Higgs boson decays into a neutralino and a chargino, where the charged Higgs 
boson is produced in association with a t ov i quark are done in [13, 20] (see also Refs. [25, 26]). 
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2 MSSM parameter space 



As noted above, Ma and tan f3 may be chosen as the MSSM inputs characterizing the MSSM 
Higgs bosons' decays into SM particles''. But when Higgs boson decays to -inos are inchidcd, 
new MSSM inputs specifying the -ino sector also become crucial. To identify the latter, the 
already mentioned Higgs/ Higgsino mixing mass, /i, and the SUSY-breaking SU{2)j^ gaugino 
mass, M2, in addition to tan/3, are required. The SUSY-breaking C/(l)y gaugino mass. Mi, 
is assumed to be determined from M2 via gaugino unification {i.e., Mi — ^tan^ 9\.yM2). This 
will fix the tree-level -ino masses (to which the radiative corrections are quite modest) along 
with their couplings to the Higgs bosons. 

Inputs (assumed to be flavor-diagonal) from the slepton sector are the left and right soft 
slepton masses for each of the three generations (selectrons, smuons, and staus) and the 
trilinear 'A-terms' which come attached to Yukawa factors and thus only Ar has a potential 
impact. A priori, all six left and right mass inputs (and A^-) arc independent. However, in 
most models currently advocated, one has ~ and ~ . Herein these equalities 
are assumed to hold. 

2.1 Experimental limits 

To maximize leptonic -ino BR enhancement, sleptons should be made as light as possible. But 
direct searches at LEP [28, 29] place significant limits on slepton masses: m~ > 99.0 GeV, 
"^^j > 91.0 GeV, m~ > 85.0 GeV (these assume that the slepton is not nearly-degenerate with 
the LSP) and m~ > 43.7 GeV (from studies at the Z'^ pole). Furthermore, the sneutrino masses 
are closely tied to the left soft mass inputs, and, to avoid extra controversial assumptions, only 
regions of the MSSM parameter space where the LSP is the lightest neutralino rather than a 
sneutrino will be considered^. To optimize the -ino leptonic BRs without running afoul of the 
LEP hmits, it is best^ to set mj^ — mj^. If all three generations have the same soft inputs 
(with At- — A^ — 0), then the slepton sector is effectively reduced to one optimal input value 



anywhere near as beneficial as those into electrons or muons, it would be even better if the 
stau inputs were significantly above those of the first two generations. This would enhance 
the -inos' BRs into electrons and muons. In the general MSSM, one is of course free to choose 
the inputs as such. Doing so would also weaken restrictions from LEP, especially for high 
values of tan p. Fig. 1 in [20] shows values for this optimal slepton mass over the M2-// plane 

''Several other MSSM inputs also enter into the radiatively-corrected MSSM Higgs boson masses and 
couplings of the MSSM Higgs bosons to SM particles, namely, inputs from the stop sector — the soft SUSY- 
breaking stop trilinear coupling At plus the stop masses — and the Higgs/Higgsino mixing mass fi. In the 
present work the stop masses are assumed to be heavy (w 1 TeV) whereas At is fixed to zero. The n parameter 
is not crucial for the SM decay modes; however, it will become so when decays to inos arc considered. 

^Further, if a sneutrino were the LSP and thus presumably the main constituent of galactic dark matter, its 
strong couplings to SM EW gauge bosons would lead to event rates probably inconsistent with those observed 
by Supcr-Kamiokande. In contrast, the coupling of an -ino to SM EW gauge bosons can be tuned to obtain 
rates consistent with current experimental limits. 

^Unless this leads to m~ < m~o < mj_^ , in which case X2 decays to charged leptons will be suppressed with 

respect to X2 decays to neutrinos, to avoid which having m~ < m~ is preferred. 



(defined 




However, since -ino decays into tau-leptons are generally not 
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relevant to the -ino sector for tan (3 = 10, 20. Setting the soft stau mass inputs 100 GeV above 
those of the other soft slepton masses, as will often be done herein, complies with current 
experimental constraints and moderately enhances the signal rates [24]. 

2.2 The signal inclusive cross sections 

Figs. 1, 2 and 3 show the LHC rates (in fb) for a{pp H^) x BR{H^ UN) + a{pp A°) 
X BR(A° — > AiN), where N is any number (including zero) of invisible neutral particles 
(in the MSSM these are either neutrinos or Xi LSPs) obtained for tan/9 = 5, 10, and 20, 
respectively^^. (Hereafter this sum of processes will be abbreviated by a{pp x 
BR{H'^, A'^ — >• UN).) Each figure gives separate scans of the ii vs. M2 plane most relevant to 
the -ino sector for (from top to bottom) Ma — 400, 500, and 600 GeV — covering the range 
of Higgs boson masses of greatest interest [24] . This is in the region of the MSSM parameter 
space where observation of h'^ alone may be insufficient to distinguish the MSSM Higgs sector 
from the SM case {i.e., the decoupling region). The darkened zones seen around the lower, 
inner corner of each plot are the regions excluded by the experimental results from LEP. 

First observe that these 'raw' or 'inclusive' {i.e., before applying selection cuts to the basic 
event-type) rates may be sufficiently large. For an integrated luminosity of 100 fb~^, the 
peak raw event number is around 4000(1700) events for Ma = 400(600) GeV and tan /5 = 20, 
irrespective of the sign of fi. Also observe that low values of and M2 yield the highest signal 
rates, though significant event numbers are also found when one but not the other of these 
parameters is increased (especially rates do fall rapidly when M2 increases much beyond 
500 GeV). These numbers are substantial (especially at high tanP) and, if experimental 
efficiencies are good, they may facilitate a much more accurate determination of some masses 
or at least mass differences in the -ino spectrum as well as the Higgs-ino mass differences than 
those achieved in previous studies based solely on H^, — > X2X2 decays. 

Note the color coding of the three figures depicting what percentage of the signal events 
are coming from Higgs boson decays to X2X2' > 90% the red zones, from 90% down to 50% 
in the yellow zones, from 50% to 10% in the blue zones, and < 10% in uncolored regions. If 
the events are not coming from X2X2^ then they are almost always from Higgs boson decays 
including heavier neutralinos, j.e^, H^,A° ^ X2X3> X2X4> XaXa^ XaX"^ X4X4 (possibly also with 
contributions from H°,A° xfx2 1X2X2 which are also taken into account here). Also note 
that the main source of events at the optimal location in the -ino parameter space shifts from 
X2X2 heavier -ino pairs as Ma grows from 400 to 600 GeV. Irrespective of the heavier Higgs 
boson masses, Higgs boson decays to X2X2 the dominant source of signal events in regions 

^"^Thcsc figiiic's arc; gcnicratcd using private codes; however, these have been cross-checked against those of 
the IS ASUS Y package of IS A JET [30] and the two are generally consistent, exceptions being a few coding 
errors in ISASUSY and the latter's inclusion of some mild radiative corrections for the slepton and -ino 
masses which are not incorporated into the codes used here. These caveats are noteworthy since results from 
the output of the ISASUSY code will be used as input for the simulation work that follows. These small 
distinctions may cause a shift in the parameter space locations of particularly-abrupt changes in the rates 
due to encountered thresholds, though the gross features found in this section and in the ISASUSY-based 
simulation studies are in very good agreement. Finally, note that higher-order corrections to the Higgs boson 
-ino -ino couplings are incorporated into neither ISASUSY nor the private code. A recent study[31] indicates 
that these generally enhance the partial decay widths by 010%; enhancement to BRs may be even more. This 
would make rates reported in this work on the conservative low side. 
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Figure 1: a{pp — > H^, A^) x BR(if°, — ^ MN) (in fb), where i = or and represents 
invisible final state particles, also showing where the percentage from H^, A^ —>■ X2X2 > 
90% (red), 50% - 90% (yellow), 10% - 50% (light blue), < 10% (white), with tan/3 = 5, 
Ma = 400 GeV (top), 500 GeV (middle), 600 GeV (bottom). Optimized slepton masses (with 
stau inputs raised 100 GeV) are used, and with rrit = 175 GeV, = 4.25 GeV, m~ = 1 TeV, 
m~ = 800 GeV, A^ = Ai = 0. The cross-hatch shaded areas are excluded by LEP. 
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Figure 2: a{pp — > ^4°) x BR(if°, — ^ MN) (in fb), -where i = or and represents 
invisible final state particles, also showing where the percentage from H^, —>■ X2X2 > 
90% (red), 50% - 90% (yellow), 10% - 50% (light blue), < 10% (white), with tan/3 = 10, 
Ma = 400 GeV (top), 500 GeV (middle), 600 GeV (bottom). Optimized slepton masses (with 
stau inputs raised 100 GeV) are used, and with rrit = 175 GeV, = 4.25 GeV, m~ = 1 TeV, 
m~ = 800 GeV, A^ = Ai = 0. The cross-hatch shaded areas are excluded by LEP. 
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Figure 3: a{pp — > ^4°) x BR(if°, — ^ MN) (in fb), -where i = or and represents 
invisible final state particles, also showing where the percentage from H^, —>■ X2X2 > 
90% (red), 50% - 90% (yellow), 10% - 50% (light blue), < 10% (white), with tan/3 = 20, 
Ma = 400 GeV (top), 500 GeV (middle), 600 GeV (bottom). Optimized slepton masses (with 
stau inputs raised 100 GeV) are used, and with rrit = 175 GeV, = 4.25 GeV, m~ = 1 TeV, 
m~ = 800 GeV, A^ = Ai = 0. The cross-hatch shaded areas are excluded by LEP. 
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with low M2 values and moderate to high values of But for low to moderate M2 values 
and low values of the dominant source of signal events shifts to the previously-neglected 
decays into the heavier -inos. Thus, inclusion of these neglected modes opens up an entirely 
new sector of the MSSM parameter space for exploration. Furthermore, the parameter space 
locations with the maximum number of signal events also shifts to these new sectors as the 
masses of the Higgs bosons rise. Therefore, the regions in MSSM parameter space wherein 
a{pp — > ^4°) X BR(i7°, ^4° — > AiN) processes can be utilized in the search for the heavier 
MSSM Higgs bosons will certainly expand substantially with inclusion of these additional 
decay channels. 

The rates illustrated in Figs. 1-3 incorporate indirect decay modes. That is, if the Higgs 
boson decays into a pair of neutralinos, and then one or both of these 'primary' neutralinos 
decay into other neutralinos (or other sparticles or the hght Higgs boson or both on- and 
off-mass-shell SM gauge bosons) which in turn give rise to leptons {with no additional colored 
daughter particles), then the contribution from such a decay chain is taken into account. This 
remains true no matter how many decays there are in the chain between the primary -ino and 
the AiN final state, the only restrictions being that each decay in the chain must be a tree- 
level decay with at most one virtual intermediate state (so 1 to 3 decay processes are included 
but not 1 to 4 decays, etc.). (As already intimated, the intermediate state is expected to be 
an on- or off- mass-shell SM gauge boson or slepton, charged or neutral.) The decay modes 
omitted due to these restrictions are never expected to be significant. Thus, effectively all 
tree-level decay chains allowable within the MSSM have been taken into account. Potential 
contributions from literally thousands of possible decay chains are evaluated and added to the 
results. 

Inspection of Figs. 1-3 supports selection of the following representative points in the 
MSSM parameter space to be employed repeatedly in this work. These are: 

Point 1. Ma = 500 GeV, tan/3 = 20, Mi = 90 GeV, M2 = 180 GeV, n = -500 GeV, 
mj = m~ = 250 GeV, m~ = m~ = 1000 GeV. 

£30 ft ^"oft ' 9 Q 

Point 2. Ma = 600 GeV tan P = 35, Mi = 100 GeV M2 = 200 GeV n = -200 GeV, 
mj - 150 GeV, m~ , = 250 GeV, m~ = 800 GeV, m~ = 1000 GeV. 

(soft ' ^so/t '9 '9 

(Also recall that nij = mj^ = mj^ and Ar = Ai = 0.) Point 1 represents a case where 

most of the signal events result from H^, A^ — > X2X2 decays^^, whereas Point 2 is a case where 
decays including heavier -inos make the dominant contribution. Here tan (3 has been set fairly 
high to enhance rates, as Figs. 1-3 suggest. 

In Fig. 4, the parameter values of Point 1 (left plot) and Point 2 (right plot) are adopted, 
save that the parameters Ma and tan/? are allowed to vary, generating plots in the Ma vs. 
tan/3 plane. Color shading on the left-side plot clearly shows that the X2X2 decay modes 
totally dominate in the production of signal events for this choice of M2, /J, -ino inputs out 
to Ma — 700 GeV. Similarly, the right-side plot shows that for the -ino inputs of Point 2 the 

^^This choice of parameters, including the degenerate soft selection, smuon and stau inputs, also corresponds 
to one of the choices adopted in [9]. 
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Figure 4: a{pp H°,A^) x BR{H°,A^ 4iN) (in fb), where £ = or jj^ and 
represents invisible final state particles for Point 1 (left side): Mi = 90 GeV, M2 = 180 GeV, 
fj, = -500 GeV, mj ^ = m~^^^ = 250 GeV, m~ = m~ = 1000 GeV; and Point 2 (right side): 

Ml = 100 GeV, M2 = 200 GeV, = -200 GeV, mj^^ J ^7,0 ft = 150/250 GeV, m~/m~ = 
800/1000 GeV. Color coding as in Figs. 1-3. 



previously neglected decay modes to heavier -inos dominate, save for a relatively small region 
around Ma ~ 350-450 GeV and tan (3 ~ 2-10. Color coding as in Figs. 1-3. 

It will be noteworthy to compare the declines in raw rates with increasing Ma and de- 
creasing tan (3 shown here to the corresponding Ma vs. tan f5 discovery region plots based on 
detailed simulation analyses presented in the analysis section to follow. 

Fig. 5 illustrates how results depend on the slepton mass(es). In the upper plot, showing 
the overall rate, a{pp — A°) x BR(iy°, A° AiN), as a function of mj ^ = mj^ ^, one 
generally sees the naively expected decline in the rate as mj increases. If the -inos decay 

^soft 

through on- or off-mass- shell sleptons, then the decay products always include leptons (and 
usually charged leptons). However, as the sleptons become heavier (first becoming kinemati- 
cally inaccessible as on-mass-shell intermediates and then growing increasingly disfavored as 
off-mass-shell intermediates), the EW gauge bosons become the dominant intermediates, in 
which case a large fraction of the time the decay products will be non-leptons, and so the BR 
to the 4i final state drops. The plot though also reveals an often far more complex dependence 
on mj , with rapid oscillations in the rate possible for modest changes in mj 

^soft ^soft 

Note again that Point 1, drawn in red in Fig. 5, represents a case where most of the signal 
events result from H^, X2X2 decays, whereas Point 2, drawn in blue, is a case where 

decays including heavier -inos make the dominant contribution. This is made clear by the 
lower plot where the percentage of the inclusive rate from X2X2 decays is plotted vs. mj 
In Fig. 5, the slepton mass is varied. But later in this work the value of m-x will be fixed at 

'^soft 
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red — Point 1 (at location of asterisk) 

blue — Point 2 (at location of asterisk) 

dashed — Point 2 but change Ma = 400 GeV, tan /3 = 30 

dotted — Point 2 but change Ma = 400 GeV, tan /? = 5 




500 



labeling same as in (a) above 
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Figure 5: Dependence on slepton mass, (a) a{pp 
where £ 



or and N represents invisible final state particles, vs 



MN) (in fb). 



f-soft 



m-r for 

''L,R 



MSSM parameter Point 1 (red) and Point 2 (blue) as well as some variations based on Point 



2 (black). Asterisks mark the mj values to be used for Points 1 and 2 later in this work. 
{b) percentage of the inclusive rate from X2X2 decays vs. m~„ 



f-soft 



with other labeling as in (a). 
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the values given earlier for Points 1. and 2. (these locations are marked by asterisks in both 
plots in Fig. 5). These choices are fairly optimal, especially for Point 1. 

Points 1. and 2. show some interesting dependence on mj . This dependence can be made 
more acute though by adjusting the input parameters. For instance, the black dotted and 
dashed curves in Fig. 5 result from lowering the value of Point 2 to 400 GeV and changing 
tan (3 from 35 to 5 and 30, respectively. Then not only docs the inclusive rate undergo rapid 
variation with mj , but the percentage of the inclusive rate from X2X2 decays fluctuates 

'-■soft 

rapidly as well. Points 1. and 2. were selected for further analysis later in this work in part 
because the results are not strongly affected by a small shift in the value of mj . However, 
apparently this is not true for all points in MSSM parameter space. 

Finally, notice that the overall normalization of both processes gg H^, and bb 
. A*^ is of 2 ^ 1 lowest-order^^. Each of these gluon- and quark- fusion partonic contributions 
is separately convoluted with an empirical set of PDFs (CTEQ 6M [34] in this case) to 
obtain predictions at the proton-proton level, for which the total center-of-mass energy is 
\/s = 14TeV. The cross-section thus defined is computed using the MSSM implementation 
[35] of the HERWIG program [36] (as available in Version 6.5 [37], with the exception of the 
choices rrit = 175 GcV and = 4.25 GeV for the top and bottom quark masses) and the 
MSSM input information produced by ISASUSY (through the ISAWIG [38] and HDECAY 
[41] interfaces). Sometimes a Higgs boson will be produced in association with jets, and thus, 
as discussed in Ref. [32], what percentage of the time a Higgs boson is produced with hadronic 
activity passing jet selection criteria (as will be applied in the analysis section) is (possibly) 
sensitive to the type of emulation (2 ^ 1 or 2 ^ 3) being employed. Note though that in 
Figs. 1-4 colored fermions are not allowed in the -ino decay chains. This is in fact inconsistent 
and leads to an over-(under-)estimate of the hadronically-quiet (inclusive, allowing jets) Ai 
rates (the under-estimation of the inclusive rates is expected to be modest due to the price of 
extra BRs in the decay chains of the neglected channels). To attempt to correct for this by 
factoring in results from the simulation runs might obscure what is meant by 'raw' rates, so 
this minor inconsistency is simply tolerated in these estimates. 

2.3 Signal-to-background rates 

The signal, taken here to be events resulting from heavy MSSM Higgs bosons decaying into 
-ino pairs, is not the only relevant quantity in this analysis that depends on the position in 
the MSSM parameter space — backgrounds from other MSSM processes will also vary from 
point to point. Fig. 1 of [27] shows the competing processes for -ino pair-production via 
Higgs boson decays^^: 'direct' -ino production {i.e., via a s-channel gauge boson) and -inos 
produced in 'cascade' decays of squarks and gluinos. The latter is considered in some detail in 

-•^^There is an alternative 2^3 approach based on MC implementation of gg/qq — > bhH^ ,bbA^ diagrams. 
The results of these two approaches have been compared and contrasted in Ref. [32] . A full MC implementation 
for the 2 — > 3 approach based on gg — > ggH^, ggA^ and related modes (eventually yielding two jets in the final 
state alongside _ff° or A'^) [33]) is as-of-yet unavailable though in public event generators. It is therefore more 
consistent to solely employ complete 2^1 emulations and not incomplete 2 —^ 3 ones. 

^^One could also consider signals from Higgs boson decays to other sparticles, especially sleptons. This 
was discussed in [39] , which demonstrated that the heavier MSSM Higgs boson decays to sleptons only have 
sufficient BRs for low values of tan/3 (~3). 
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[27], but will be removed from consideration here by making the assumption throughout this 
work that gluinos and squarks are heavy {circa 1 TeV). However, since the signal depends on 
them, (all) the -inos cannot be made heavy^^, and the masses of the EW gauge bosons are 
known, so the direct channel background cannot be easily removed by restricting the analysis 
to some subset of the parameter space by means of such a straight-forward assumption. 

In fact, the location in the parameter space where the raw signal rate is largest sometimes 
differs from that where the ratio of the signal to the leading background from direct -ino 
production is largest. For instance, the plot in Fig. 2 (tan/3 = 10) for Ma — 600 GeV shows 
a maximum in the inclusive rate at approximately {/j,M2) = (—200 GeV, 250 GeV). On 
the other hand, the signal-to-background ratio {S/B) is largest at (—250 GeV, 500 GcV). 
The production cross-section for the Higgs bosons is the same at both points. Thus, to 
understand why the two locations differ so much the BR(i7°, — > MN) and the direct -ino 
production x BR( -inos — > MN) need to be studied. The former drops from ~6% to ~2% 
in moving from the inclusive rate maximum to the S/B maximum (thus cutting the overall 
signal rate by a factor of 3). The background at the inclusive rate maximum is mostly X2X3) 
X3X4 and xt)(2 with respective production cross-sections (and BRs into MN final states) of 
4 X 10-2 pb (18%), 1 X 10-2 pb (8%) and 1 x lO-^pb (2%). At the point where the S/B is a 
maximum, these (still dominant) backgrounds rates shift to 1 x lO-^pb (16%), 1 x 10-^ pb 
(27%) and 1 x lO-^pb (2%), respectively. So the X2X3 production rate drops by a factor 
of 4 while XsXd production almost vanishes (which is the main factor), mostly because of 
increased phase space suppression due to larger -ino masses: m~o(m~o)[m~o]{m~±} changes 
from 118(180)[212]{289}GeV at the rate maximum to 219(257)[273]{515} GeV at the S/B 
maximum. The result is that the overall background rate drops by a factor of 5. In short, the 
S/B improves because the direct -ino pair-production cross-section falls more rapidly than 
the signal BR into AiN final states. Analogous plots to those in Figs. 1-3 studying the S/ B 
variation across the parameter space are not presented. Instead, discovery regions for selected 
-ino input parameter sets will be given in Sect. 4. While favorable MSSM points have been 
chosen for the simulation analyses, they were not selected to maximize the S/B. Therefore, 
this channel may work even better at points other than those analysed in detail herein. 

3 mSUGRA parameter space 

Augmenting the general MSSM with additional assumptions about the unification of SUSY 
inputs at a very high mass scale yields the more restrictive 'mSUGRA' models. Here the 
number of free input parameters is much reduced (hence the popularity of such scenarios 
for phenomenological analyses), with said free parameters generally set as tan/3, a universal 
gaugino mass defined at the Grand Unification Theory (GUT) scale (M1/2), a universal GUT- 
level scalar mass (Mq), a universal GUT-lcvcl trilinear scalar mass term {Aq), and the sign of 
fi (henceforth, sgn(/i)). As ah^cady noted, the signal has a strong preference for low values of 
Yet in mSUGRA scenarios, is not a free parameter, as it is closely tied to the masses of 
the scalar Higgs bosons via the Mq input. An earher study of charged Higgs boson decays into 

^^The sleptons also cannot be made arbitrarily heavy. Direct slepton pair-production, as studied in [40], 
will generally lead to dilepton final states rather than the 41 final state desired here. The smaller contributions 
from these processes are included in the analyses to follow. 
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Figure 6: a{pp H^,A^) x BR(if°,A° 4iN) (in fb), where i = or fi^ and 
represents invisible final state particles for tan/5 = 5,10,20 in the mSUGRA Mq vs. M1/2 
plane, with sgn(yu) = +1 and Aq = 0. Colors depict the percentage of events stemming from 
A° X2X2 > 90% (red), 50% - 90% (yellow), 10% - 50% (light blue), < 10% (white). The 
dark shaded regions are excluded by theoretical considerations or LEP measurements (save 
constraints from LEP Higgs-strahlung which roughly reach up to the dashed green curves with 
considerable uncertainty — see text). Also shown in purple are the CMS TDR (Fig. 11.32) 5a 
discovery regions (assuming L^nt = 30 fb^^) for H^,A^ X2X2- The solid purple lines show 
the extent of the plots in Fig. 11.32. 
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a neutralino and a chargino [20] demonstrated that this was sufficient to preclude detection of 
a,3£ + top-quark signal from such processes over the entire reach of the unexcluded mSUGRA 
parameter space. Here, with the heavier neutral MSSM Higgs bosons, the situation is not so 
discouraging. Fig. 6 shows the values for a(j)p — > H^,A^) x BR(if°, A'^ — > MN) obtained for 
tan (3 = 5, 10, 20 and fi > 0. Two disconnected regions of unexcluded parameter space appear 
where the expected number of events (for 100 fb~^ of integrated luminosity) is in the tens to 
hundreds (or even thousands). Interestingly, one of these (which includes discovery regions 
depicted in the CSM TDR [14]^^) is where X2X2 is the dominant source of M events while the 
other is where decays of the heavier -inos dominate. For tan [3 = 5, rates in the X2X2 region 
are much larger than in the heavier -inos region. However, for tan/3 = 20, rates in the two 
regions become more comparable. 

Also shown as sohd purple zones on the tan (3 = 5 and tan /3 = 10 plots are 5a discovery 
regions from the CMS TDR (Fig. 11.32) [14]. These CMS TDR discovery regions assume an 
integrated luminosity of just 30 fb""*^, and thus would have certainly been considerably larger 
if a base luminosity of 100 fl:)^^ was used instead. This CMS TDR analysis was at a technical 
level comparable to that in this work, hut only considered MSSM Higgs boson decays into X2X2 
pairs. Thus, the CMS TDR analysis would not pick up the region where heavier -ino decays 
dominate (in fact the plots in Fig. 11.32 in the CMS TDR only showed the regions dehneated 
by the solid purple lines in Fig. 6). Given that the somewhat lower rates of the higher Mq, 
heavier -ino decays-dominated region may be compensated by assuming a larger integrated 
luminosity, as well as perhaps finding a higher selection efficiency due to harder daughter 
leptons, it is difficult to infer from the CMS TDR 5a 30 fb~^ discovery regions whether or not 
disjoint discovery regions may develop in this novel region of the parameter space. This is 
currently under investigation [15]. 

The excluded regions shown in Fig. 6 merit some explanation. Note that in each plot the 
discovery region from the CMS TDR cuts into the excluded region, whereas in Fig. 11.32 
of the CMS TDR they do not touch the (more limited) excluded regions shown. This is 
mainly because the excluded regions in Fig. 11.32 of the CMS TDR only mark off regions 
where the Xi is not the LSP (because the mass of the lighter stau is lower — this removes 
the upper left corner of the plots) and where EWSB is not obtained (along the horizontal 
axis), while ignoring other experimental constaints — such as the lower limit on the lighter 
chargino's mass from the LEP experiments. Such additional experimental constraints are 
included, for instance, in the excluded regions shown in Fig. 20-1 of the ATLAS TDR [16] 
These experimental constraints have been updated to represent the final limits from the 
LEP experiments, accounting for the gross differences between the excluded regions depicted 

^^Note: virtually all mSUGRA parameter space plots in the TDR showing excluded regions are for tan/3 ~ 
10; the exceptions being the tan/3 = 5 plot in Fig. 11.32 and the tan/3 = 35 plots in Figs. 13.12 & 13.13; 
and the tan/3 = 35 plots seem to inaccurately have the tan/3 = 10 exclusion zones. These latter plots and 
others in Chapter 13 do show a chargino lower mass limit (green dotdashed curve) and other supercollider 
experimental bounds which are more consistent with the excluded regions shown in the ATLAS TDR (and in 
the present work). 

^^Notc: virtually all mSUGRA parameter space plots in the TDR showing excluded regions are for tan /3 = 10 
and for (the now ruled-out) tan/3 = 2. 
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in the ATLAS TDR and those in the present work^^. Somewhat crude^^ estimates for the 
regions excluded by the LEP searches for MSSM Higgs bosons are indicated separately by the 
dashed green hnes based on the empirical formula developed by Djouadi, Drees and Kneur [18]. 
Finally, it must be emphasized that constraints from lower-energy experiments (in particular 
from b s'-f) and from cosmological considerations (such as LSP dark matter annihilation 
rates) are not herein considered. In the far more restricted parameter domain of mSUGRA 
models it is more difficult to circumvent such constraints, and they can exclude considerable 
portions of the allowed parameter space shown in the figures (for further details, see [19]). 

As was done with the general MSSM parameter space, Fig. 5 enables selection of a couple 
of representative mSUGRA points for simulation studies. These are: 

Point A. Mo = 125 GeV, M1/2 = 165 GeV, tan/3 = 20, sgn(/i) = +1, Aq = 0. 
Point B. Mo = 400 GeV, M1/2 = 165 GeV, tan/3 = 20, sgn(/i) = +1, Aq = 0. 

Point A is dominated by A^ — » X2X2 ~^ 4£ decays (which account for more than 99% of 
the inclusive signal event rate before cuts) while in Point B the corresponding rates are below 
30% (the largest signal event channel is now H^,A° X1X2 4£, yielding over 50% of the 
events, with significant contributions from H^,A^ — > X2X2) X2X3) X2X4 ~^ 

U). Full MC and 

detector simulations for Points A and B will be presented in the next section. These will show 
that AIN signals remain visible in the mSUGRA parameter space, at least at these points. 

4 Simulation analyses 

The HERWIG 6.5 [37] MC package (which obtains its MSSM input information from ISASUSY 
[30] through the ISAWIG [38] and HDECAY [41] interfaces) is employed coupled with private 
programs simulating a typical LHC detector environment (these codes have been checked 
against results in the literature). The CTEQ 6M [34] set of PDFs is used and top and bottom 
quark masses arc set to uit = 175 GeV and = 4.25 GeV, respectively. 
Four-lcpton events are first selected according to these criteria: 

• Events have exactly four leptons, £ — e or /i, irrespective of their individual charges, 
meeting the following criteria: 

""^^Raising the lower bound on the chargino mass from the circa 1998 [17] LEP-1.5-era ^65 GeV to ~100 GeV 
raises the approximately horizontal boundary for higher Mq values, while the rise of the bounds for the slepton 
masses from -^45 GeV to m~ , m~ , m~ ~ 99 GeV, 91 GeV, 85 GeV adds the quarter-circle-like bite seen in the 
lower-left corner of the tan/3 = 10 plot in Fig. 6 (which is absent in the ATLAS TDR plots). 

^^For reasons detailed in [18], foremost among which is the uncertainty in the calculation of M/,. Herein the 
Higgs boson mass formulas of ISAJET [30] and [12] are employed. Results here are roughly consistent with 
Figs. 1 & 2 of [18] (2006 paper). Note that in the case of mSUGRA, unlike in the general MSSM examples in 
the current work, the stop and other squark parameters — which make the main contributions to the quite 
significant radiative corrections to Mh — are determined from the few mSUGRA inputs without the need to 
set values by hand for assorted soft SUSY-brcaking masses. Certainly, in mSUGRA, the LEP bounds on light 
Higgs boson production are strongly-tied to rates for heavy Higgs boson to sparticle decay channels, though 
this correlation will not be intensively examined in this work. 
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Each lepton must have \ri^\ < 2.4 and E^>7,4 GeV for e,/x (see ATLAS TDR [16]). 
Each lepton must be isolated. The isolation criterion demands there be no tracks (of 
charged particles) with px > 1.5 GeV in a cone of r = 0.3 radians around a specific 
lepton, and also that the energy deposited in the electromagnetic calorimeter be less 
than 3 GeV for 0.05 radians < r < 0.3 radians. 

Aside from the isolation demands, no restrictions are placed at this stage on the amount 
of hadronic activity or the number of reconstructed jets in an event. 

Further, 

• Events must consist of two opposite-sign, same-flavor lepton pairs. 

Events thus identified as candidate signal events are then subjected to the following cuts: 

• Z°-veto: no opposite-charge same-flavor lepton pairs may reconstruct Mz ± 10 GeV. 

• restrict Ej,: all leptons must flnally have 20 GeV < Ej, < 80 GeV. 

• restrict missing transverse energy, i?™'®^: events must have 20 GeV < E'™**^ < 130 GeV. 

• cap all jets must have < 50 GeV. 

Jets are reconstructed using a UAl-like iterative {i.e., with splitting and merging, see 
Ref. [42] for a description of the procedure) cone algorithm with fixed size 0.5, wherein 
charged tracks are collected at £'r > 1 GeV and \r]\ < 2.4 and each reconstructed jet is 
required to have E^^^ > 20 GeV. 

Lastly, application of an additional cut on the four-lepton invariant mass is investigated: 

• four-lepton invariant mass (inv. m.) cut: the 4i inv. m. must be < 240 GeV . 

For the signal events, the upper limit for the four-lepton inv. m. will be Mh,a — 2M~(), and 
thus its value is dependent upon the chosen point in MSSM parameter space. In the actual 
experiment, the value of Mh a — 2M~o would be a priori unknown. So one could ask how a 
numerical value can be chosen for this cut? If too low a value is selected, many signal events 
will be lost. On the other hand, if too large a value is chosen, more events from background 
processes will be accepted, diluting the signal. One could envision trying an assortment of 
numerical values for the four-lepton inv. m. upper limit (one of which could for instance be the 
nominal value of 240 GcV noted above) to see which value optimized the signal relative to the 
backgrounds. However, here sparticle production processes are very significant backgrounds 
(after application of the other three cuts, only such processes and residual events 
remain), which, like the signal, may well have unknown rates. Thus, strengthening this cut 
would lower the total number of events without indicating whether the signal to background 
ratio is going up or down — unless additional information is available from other studies 
at least somewhat restricting the location in MSSM parameter space Nature has chosen. If 
such information were available, this cut could indeed lead to a purer set of signal events. 
One could instead consider all events from MSSM processes to be the signal while the SM 
processes comprise the background. However, the aim of this work is to identify the heavier 
Higgs bosons, not merely to identify an excess attributable to SUSY. 
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Table 1: Relevant sparticle masses (in GeV) for specific MSSM and mSUGRA parameter 
points studied in the analyses. 





Point 1 


Point 2 


Point A 


Point B 


Ma 


500.0 


600.0 


257.6 


434.9 


Mh 


500.7 


600.8 


257.8 


435.3 


Xl 


89.7 


93.9 


60.4 


60.8 


~i\ 
X2 


176.3 


155.6 


107.8 


108.0 


— n 

xl 


506.9 


211.8 


237.6 


232.8 


xi 


510.9 


262.2 


260.0 


256.3 


xt 


176.5 


153.5 


106.8 


106.8 


xt 


513.9 


263.2 


260.0 


258.2 


m~ 


241.6 


135.5 


154.8 


407.9 




253.8 


156.3 


145.7 


406.1 


^Hi 


252.0 


154.3 


145.6 


406.1 


^e2 


254.4 


157.2 


174.1 


415.7 




256.2 


159.2 


174.2 


415.7 


m~ - TO- 


0.59 


0.96 


28.46 


9.56 




4.20 


4.81 


28.62 


9.63 



Detailed results are tabulated for the aforementioned two general MSSM and two mSUGRA 
parameter space points. MSSM Point 1 and mSUGRA Point A have the vast majority of their 
4£ events from A° X2X2, while MSSM Point 2 and mSUGRA Point B obtain most of 
their 4i events from Higgs boson decays to heavier -ino pairs (X2X3) X2X4 X3X3) X3X4 and/or 
X4X4)- The sparticle spectra^^ for these points are presented in Table 1. 

4.1 MSSM benchmark points 

Table 2 shows results for MSSM Point 1, a H^, X2X2'^'^^^^^^^^ point. Note that, after 

cuts, signal events do make up the majority of events in the sample. The only remaining 
backgrounds are from direct neutralino/chargino pair-production^° (denoted by xx), from 
slepton pair-production (denoted by i, V) and from production. The number of 

events obtained from decays after cuts is about twice the number obtained from decays. 
This is despite the fact that the and A^ production cross sections are the same within 
1%. The ratio of A^ to events at this point can be compared to that for inclusive rates 

-•^^The older ISASUSY version which inputs sparticle masses into HERWIG 6.3 lacks D-terms in the slepton 
masses, meaning the smuon masses in the simulation runs equate to the selectron masses given in Table 1. 
This has a minor effect upon the edges in the Dalitz-like 'wedgebox' plots to be shown later. See discussion 
in [27]. 

■^^Herein final states involving a sparton and a chargino/neutralino are included together with the results 
for XX) as designed in HERWIG, though for the points studied here the latter overwhelmingly dominate the 
former. 
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Tabic 2: Event rates after the successive cuts defined in the text for MSSM Point 1 (assuming 
an integrated luminosity of 100 fb~^). 



Process 


M events 


£+£-£(0+^(0- 


Z°-veto 




TTimiss 


Ef 


M inv. m. 


g , 9 


118 


64 


49 


19 


1 








e,v 


100 


65 


46 


30 


23 


13 


7 


XX, (l/9X 


34 


17 


13 


10 


5 


2 


1 


tH' + c.c. 

























1733 


1683 


43 


39 


5 


4 


4 




47 


23 


2 


1 


1 










1 


2 


2 


1 


1 








A° signal 


20,32 


18,31 


14,26 


13,25 


11,22 


8,17 


6,13 



(with no cuts) which may be calculated using the BRs obtained from ISASUSY^^. Including 
all possible decay chains, ISASUSY numbers predict : = 1.83 : 1.00 (64.7% events). 
This is in reasonable agreement with A^ : — 1.6 : 1.0 (61.5% A^ events) obtained from 
the 4£ before cuts entries in the first column of Tabic 2. The different and ^4° event 
rates may then be traced back to differences in the / A^-X2-X2 couplings (as opposed to 
the enhancing or opening up of other decay modes, such as for instance h^h^). 
Study of the inclusive rates based on the ISASUSY BRs also confirmed that over 99% of the 
four-lepton signal events resulted from H^/A^ X2X2 decays. The percentage of A^ — > A£ 
events surviving the subsequent cuts is about 10% larger than the percentage of — > M 
events surviving. 

Fixing the -ino input parameters M2 & /i and the slepton & squark inputs to be those of 
MSSM Point 1, tan/? and Mj^ were then varied to map out a Higgs boson discovery region in 
the traditional (M^, tan (3) plane. This is shown in red in Fig. 7, where the sohd (dashed) red 
border dehneates the discovery region assuming an integrated luminosity of 300 fb~^ (100 fb~^). 
The exact criteria used for demarcating the discovery region is that there be at least 10 signal 
events and that the 99%-confidence-level upper limit on the background is smaller than the 
99%-confidence-level lower limit on the signal plus background. Mathematically, the latter 
condition translates into the formula [43]: 



^signal > (2-32)^ 



, ^V^bckgrd 

1 + — 

2.32 



(2) 



where A'gjgj^g^j and A/^j^^^i^gj-f^ arc the expected number of signal and background events, re- 
spectively. As with MSSM Point 1, direct ncutralino/chargino pair-production, slepton pair 
production and SM Z'^'^*^Z^^*^ are the only background processes remaining after cuts (the ac- 
tual number of surviving background events varies modestly with tan (3) at all points tested. 



These were normalized using HERWIG production cross-sections, though here this is of scant importance 
since the and A" production cross-sections are almost the same. Also, for consistency with the HERWIG 
simulation analysis, ISASUSY Version 7.56 was used to generate the BRs. 
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with slepton pair production continuing as the dominant background. Taking into account 
these backgrounds, 24-28 (38-45) signal events are required to meet the criteria for 100 fb~^ 
(300 fb~^) of integrated luminosity, depending on the value of tan/3, if the four-lepton inv. 
m. cut is not employed. Adding in this last optional cut changes the required numbers to 
19-22 (28-34) signal events and shifts the discovery region boundaries to those shown as blue 
(dashed blue) curves in Fig. 7. This places MSSM Point 1 just outside the upper Ma edge of 
the 100 fb~^ discovery region (whether or not the four-lepton inv. m. cut is used). Lowering 
Ma to 400 GeV raises the number of signal events from 25 to 36. Note that Fig. 4 (left-side 
plot) predicts that i?", A" decays will generate the bulk of the signal throughout the 

discovery region. The lower Ma edge of the discovery region closely follows where the (dom- 
inant) X2X2 decay becomes kinematically accessible, z.e.. Ma > 2m~o. The contribution 
outweighing the contribution was found to be a general result valid for almost^^ all points 
in the {Ma, tan /3)-plane tested: events from equaled or outnumbered those from H^. Note 
from Table 2 that MSSM Point 1 at = 500 GeV and tan /3 = 20 yielded A° : //° = 2.1 : 1.0 
(68% A^ events) after all cuts save the four-lepton inv. m. cut (as comparison to the numbers 
in the preceding paragraph indicate, A^ events tend to do slightly better at surviving the cuts, 
though little reason could be found for this small effect). Lowering Ma to 400 GeV shifts this 
ratio to A^ : ^ 3.9 : 1.0 (81% A° events). 

The preponderance of A^ events is generally greatest for lower values of Ma- For 
Ma < 375 GeV, 90-100% of the signal events are from Since Ma < Mh and Ma ^ 2m~o 
this is mainly a threshold effect. The A^ event percentage drops to around 70% when 
Ma — 415 GeV. For higher Ma values inside the 100 fb^ discovery region (outside the 
100 fb~^ discovery region but inside the 300 fb""*^ discovery region), this percentage ranges 
from ~70% down to ~55% (~60% down to ~50%), save for the upper tip where tan/3>30 
wherein the A^ percentage remains above 70% or even 80%. 

Inclusion of the four-lepton inv. m. cut with the nominal cut-off value of 240 GeV shifts 
the discovery region boundaries in Fig. 7 from the red curves to the blue ones. There are 
slight gains for low Ma values at high and low values for tan/?; however, the high Ma edges 
also recede somewhat. Note also that the highest and lowest tan /3 values which fall inside 
the discovery region are virtually unaltered. Though the cut's effect on the expanse of the 
discovery region is quite modest, inclusion of this cut at included points with lower Af^ values 
can certainly raise the signal : background. For instance, at (M^, tan/5) = (400 GeV, 20), 
this ratio goes from 37 : 19 without the 4£ inv. m. cut to 37 : 12 with it. However, shifting 
Ma to 500 GeV as in MSSM Point 1 is enough to remove any advantage, as can be seen in 
Table 2. 

Input parameters for MSSM Point 1 were also chosen to match a point studied in a previous 
analysis [9] — which only looked at X2X2 Higgs boson decays^^. The light purple contour shown 

^^Insidc of the discovery region (for 300 fb~^), a couple points along the high Ma - lower tan/J edge were 
found where the rate from very slightly exceeded that from . 

A different simulation of the quark- fusion channel involving b (anti)quarks (in the CMS note the simulation 
was performed using gg bbH^ .hhA'^ ) is adopted here. In addition, the MC analysis in [9] was done with 
PYTHIA version 5.7 [44], which only implemented an approximated treatment of the SUSY sector, while herein 
ISASUSY is used in conjunction with HERWIG (though intrinsic differences between the two generators in the 
implementation of the PS and hadronization stages should be minimal in our context). Also, the background 
processes tH~ + c.c, tiZ and tt/i, which were not emulated in [9], in this study were checked to yield no 
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Figure 7: Discovery region in red in (M^,tan/?) plane for -ino/slepton parameters ^ = 
-500 GeV, M2 = 180 GeV, Mi = 90 GeV, mj = m~ = 250 GeV as in MSSM Point 
1 (whose location is marked by a black asterisk). Here Higgs boson decays to X2X2 totally 
dominate. Solid (dashed) red border delineates the discovery region for Li^t = 300 fb~^ 
(100 fb"^). The two green curves are Ma, MH—2m~o. Also shown in light purple are analogous 

results from a previous study [9] for 100 fb^^. The blue contours add the extra cut on the 
four-lepton inv. m. for the nominal cut-off value of 240 GeV. 



background events throughout Fig. 7. 
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Table 3: Percentage of H^, 4iN events (excluding cuts) coming from various -ino 

channels for MSSM Point 2. (Other channels are negligible.) 



- xlxl 31.5% 
A' - xlxl 31.1% 
A' ^ xtxl 13.4% 

- xlxl 8.4% 
^ xf)S 6.9% 
^ 4.3% 
^ Xlxl 1.9% 

- X^X^ 0.8% 
^ X^X^ 0.75% 
^ X^X2" 0.6% 



all other contributions < 0.5% 



in the plot is the result from this older study (see the blue contour in Fig. 19 therein). Results 
in the present case for the most part agree with those of that previous study, though in the 
current analysis the discovery region extends to somewhat higher values of and dies for 
tan/3 values below ~5. The latter is primarily due to low tan/3 strong enhancement of the 
H^(A'^)-t-t coupling, which is proportional to esc p {cot j3), increasing the H^, ^4° — >■ tt BRs at 
the expense of the -ino BRs^^ BR(i/° tt) (BR(A° tt)) rises from around 0.30 to 0.68 
to 0.93 (0.51 to 0.79 to 0.96) as tan/3 runs from 6 to 4 to 2. 

For MSSM Point 2, Higgs boson decays to the heavier neutralinos and charginos neglected in 
previous studies produce most of the signal events. Table 3 gives the percentage contributions 
to the signal events among the , decay modes based on an inclusive rate study using 
BR results from ISAJET (ISASUSY) 7.58 normalized with HERWIG cross-sections. This 
parton-level analysis merely demands exactly four leptons in the (parton-level) final state. 
According to this inclusive rates study, Higgs boson decays to X2X2 contribute less than 
one hundredth of one percent of the signal events, in stark contrast to MSSM Point 1 where 
such decays accounted for virtually all of the signal events. Applying all the cuts at the full 
event-generator level does not alter this. Said numerical results with the application of the 
successive cuts for MSSM Point 2 are given in Table 4.1. 

Note that the four-lepton inv. m. cut, with the nominal numerical value of 240 GeV, removes 
about 74% of the signal events while only slightly reducing the number of background events. 
This clearly shows that this cut, while helpful for points with lower Ma values in Fig. 7, is 
quite deleterious at MSSM Point 2. Without the 4£ inv. m. cut, an integrated luminosity of 
25fb~^ is sufficient to meet the discovery criteria; while with the 4£ inv. m. cut, an integrated 
luminosity of ~130 fb^^ is required. Choosing a higher numerical cut-off would lead to a viable 

^''Thc partial widths for and decays to -inos also drop by roughly a factor of 2 in going from tan /3 = 6 
to tan/j = 2 (at Ma = 450 GeV), and the h?h^ and h'^Z'^^*^ widths increase by about a factor of 

2. These also lower the signal rate. On the other hand, decay widths to /;-quarks and tau-leptons also drop 
by a bit over a factor of 2, helping the signal. These effects are overwhelmed by an almost order-of-magnitude 
enhancement in the and A° to tt decay widths. 
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Table 4: Event rates after successive cuts as defined in the text for MSSM Point 2 
(assuming 100 fb~^). 



Process 


M events 


£+£-£(')+£(')- 


Z°-veto 




pmiss 

Jjjrp 




A£ inv. m. 


9 


817 


332 


197 


96 


21 








e,v 


12 


5 


4 


4 


2 


2 


2 


XX, (i/gx 


123 


74 


32 


17 


13 


10 


4 


tH~ + c.c. 


76 


38 


22 


15 


9 


3 


1 




1733 


1683 


43 


39 


5 


4 


4 


ttz°* 


47 


23 


2 


1 


1 








tth^ 


4 


1 


1 


1 


1 








H^, AO signal 


189,179 


156,149 


64,80 


55,64 


43,50 


32,37 


9,9 



cut for this point; however, it may prove impossible to a priori decide on an appropriate value 
for the actual experimental analysis (see earlier discussion). 

Table 4.1 gives a ratio of 4£ events to H'^ A£ events (before additional cuts) of 

A° : = 1 : 1.05 (48.6% events). ISASUSY BR studies of the inclusive four-lepton event 
rates at this point also predict that will produce more signal events than ^4° this time, with 
^4° : = 1 : 1.36 (42.4% A^ events). Exact agreement between the two methods is certainly 
not expected, and it is at least reassuring that both predict more 4i events (unlike at 

MSSM Point 1). The percentage of A^ 4i events surviving the subsequent cuts is again 
slightly larger than that for H'^ — > 4i events (21% vs. 17%, excluding the four-lepton inv. m. 
cut) . Note that the Z°-veto takes a larger portion out of the signal event number for MSSM 
Point 2 than it did for MSSM Point 1, with only about 50% surviving for the former while 
about 80% survive for the latter. This is understandable since, for MSSM Point 1, virtually all 
events were from X2X2 pairs, and X2 is not heavy enough to decay to Xi via an on-mass-shell 
For MSSM Point 2, on the other hand, a variety of heavier -inos are involved, and the 
mass differences between Xs or Xa and Xi do exceed Mz- 

Again the -ino input parameters M2 & /i and the slepton &: squark inputs are fixed, this 
time to be those of MSSM Point 2, and tan/? and M4 allowed to vary to map out the Higgs 
boson discovery region in the {Ma, tan/?) plane (using the same criteria as in Fig. 7) shown 
in red in Fig. 8. As before, the solid (dashed) red border delineates the discovery region 
assuming an integrated luminosity of 300 fb~^ (100 fb~^). Assuming that the four-lepton inv. 
m. cut is omitted, MSSM Point 2 hes firmly inside the 100 fb~^ discovery region (with the 15 
sparticle/charged Higgs boson + 4 event background. Relation (2) requires 26 signal 

events to be included in the 100 fb^^ discovery region, while 69 signal events are expected). 
Note that Fig. 4 (right-side plot) predicts that — > X2X2 decays will only generate a 
substantial number of signal events when tan /3 and Ma are small (the red and yellow zones in 
the plot), with decays to heavier -inos dominating elsewhere. This leads to a disjoint discovery 
region in Fig. 8, consisting of a smaller mainly X2X2"dominated portion for lower values of 
tan f3 and Ma and a novel larger portion at considerably higher Ma values that stretches up to 
tan (3 values well above 50. Note the distance between the lower Ma edge of this larger portion 
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of the discovery region and the curves for Ma, Mh — 2m~o. In concurrence with the percentage 
contributions for MSSM Point 2 given above, the lower Ma edge of the discovery region abuts 
the Ma, Mh — m~o — m~o curves (shown in green in Fig. 8), /or tan /3> 10. The situation for 
tan/3<10 and 450 GeV < Myi <700 GeV (in both the upper and lower disjoint portions of the 
discovery region) is more complicated, with X2X2 several other decays making significant 
contributions. 

The discovery region shown in Fig. 8 represents a significant extension of LHC MSSM Higgs 
boson detection capabilities to quite high Higgs boson masses. With 300 fb"^ of integrated 
himinosity, there is some stretch of Ma values covered for almost all values of tan/? (1 < 
tan/? < 50), the exception being 4<tan/?<6. If the integrated luminosity is dropped to 
100 fb~^, the higher Ma portion of the discovery region recedes up to tan /?> 8-10, still lower 
than the 300 fb~^ discovery regions from MSSM Higgs boson decays to third generation SM 
fermions found in the ATLAS [45] and other [46] simulations. The new discovery region has 
considerable overlap with the so-called decoupling zone, where the light MSSM Higgs boson 
is difficult to distinguish from the Higgs boson of the SM, and, up to now, no signals of the 
other MSSM Higgs bosons were known. 

Though the number of signal events swells to over 50 (30) per 100 fb"^ for tan/3<2 (4), the 
background from -ino pair-production via EW gauge bosons is also becoming quite large, and 
thus more integrated luminosity is required for the excess from Higgs boson decays to meet 
the (2) criterion. Note how an 'excess' attributed to the Higgs boson signal could alternatively 
be accounted for by the MSSM background if the value of tan f3 is lowered. (Note also though 
that restrictions from LEP experiments exclude the most sensitive region of extremely low 
tan/3 values.) As in Fig. 7, the low Ma edge of the lower portion of the discovery region in 
Fig. 8 abuts the Ma, Mh — 2m~o curves. 

Yet for Ma in the vicinity of 350 GeV to 450 GeV, the discovery regions in Fig. 7 and Fig. 
8 resemble mirror images of each other: the former lies exclusively above tan/? ^ 5 while 
the latter lies exclusively below tan/3 ~ 5. The reasons behind this stark contrast, though 
a bit comphcated, critically depend on the different inputs to the slepton sector. In Fig. 8, 
for Ma ~ 470 GeV, Higgs boson decays to other heavier -inos are kinematically inaccessible, 
and, for higher tan/3 values, X2 decays almost exclusively via sneutrinos into neutrinos and 
the LSP, yielding no charged leptons. This is not the case in this region of Fig. 7 — here 
X2 undergoes three-body decays via off-mass-shell sleptons and Z^* with substantial BRs 
into charged leptons. The situation for Fig. 8 changes as tan /3 declines below ~10 since X2 
BRs to charged sleptons, while still much smaller than those to sneutrinos, grow beyond the 
percent level — sufficient to generate a low tan/3 discovery region in Fig. 8. One might expect 
analogous behavior in Fig. 7; however, in the low tan/3 region of Fig. 7 the partial widths 
r{H^,A'^ — > X2X2) much smaller, especially for A^, than they are in this region of Fig. 8 
and dechne with falling tan /3, whereas in Fig. 8 T{H^ — > X2X2) actually increases (though only 
moderately) as tan /3 falls. The X2X2 partial widths coupled with the subsequent X2 decays to 
charged leptons are large enough in the case of Fig. 8 so that the signal is not overwhelmed 
by the rising r{H^,A^ — > tt) partial widths as it is in the case of Fig. 7. Also, in Fig. 8 but 
not in Fig. 7, as Ma increases beyond ~450 GeV, contributions from other -ino pairs besides 
X2X2 become significant and further enhance the low tan /3 4£ signal rate. 

Differences in the discovery regions at very high tan /3 values are also attributable to the 



23 



CO. 




300 



400 



500 



600 



700 



800 



900 



M.(GeV) 



Figure 8: Discovery region in red in (Myi,tan/?) plane for -ino/slepton parameters /i = 
-200 GeV, M2 = 200 GeV, Mi = 100 GeV, mj = 150 GeV, m~ , = 250 GeV as in MSSM 

^soft soft 

Point 2 (whose location is marked by an black asterisk). Here Higgs boson decays to a variety 
of higher mass -inos (see text) constitute the majority of the signal events. Solid (dashed) 
red border delineates the discovery region for Lint = 300 fb"^ (100 fb"^). The green curves are 



Ma,Mh 



m-om-o and Ma, 



2,3,4; k = 1,2). The blue contours add 



the extra cut on the four-lepton inv. m. for the nominal cut-off value of 240 GeV. 



slepton input parameters. In Fig. 8, the discovery region reaches up well beyond tan/? = 50, 
while in Fig. 7 the discovery region is curtailed, ending before reaching tan f3 = 35. Since the 
soft slepton mass inputs for all three generations are degenerate for MSSM Point 1, for high 
tan P values in Fig. 7 splitting effects with the staus drive one of the physical stau masses 
well below the selectron and smuon masses. This leads to lots of -ino decays including tau 
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leptons, virtually shutting down the decays to electrons and muons. Since the soft stau mass 
inputs are elevated well above the other slepton inputs for MSSM Point 2, this high tan /? cap 
is removed in Fig. 8. 

Comments made above for MSSM Point 2 about the increased severity of the Z'^-linc cut 
and the inappropriateness of the four-lepton inv. m. cut (with the numerical cut-off set to 
240 GeV) are also applicable to points throughout the larger portion of the discovery region of 
Fig. 8. As can be seen from the blue curves in Fig. 8, inclusion of the 240 GeV M inv. m. cut 
ehminates about half of the 300 fb~^ discovery region and far more than half of the 100 fb~^ 
region, including all points between tan /3 ~ 8 and tan /3 ~ 25 for the latter. 

Also, in contrast to the discovery region of Fig. 7, in large segments of the Fig. 8 discovery 
region the number of signal events from decays exceed those from A'^ decays. First consider 
the smaller, low tan/?, portion of the disjoint discovery region. Herein, to the right of the 
Ma, Mff—m-o —^^o curves (shown in green in Fig. 8), the percentage of ^4° events ranges from 
~30-~40% (~25-~30%) for tan/3 > 2 (< 2). To the left of these curves, the A" event percentage 
grows to ~45-~60% for tan f3>2; increasing further to ~70-~80% near the region's upper left 
tip {Ma in the vicinity of 350 GeV and tan/3 around 3 to 4.5) where the signal is dominated 
by Higgs-mediated X2X2 production. 

In the novel and larger high tan/3 portion of the discovery region in Fig. 8, where the 
X2X2 contribution is minor to insignificant, the and ^4° contributions to the signal events 
stay within 20% of each other (with the A° event percentage ranging from ~40-~60%) to 
the right of the Ma, Mh — 2m~o curves. In the finger-like projection between the nearly- 
vertical Ma, Mh — m~o — m~o and Ma, Mh — 2m~o curves the A° percentage drops to < 25% 

A3 A,4 A4 

(after cuts, excluding the 4£ inv. m. cut)^^, meaning that the number of events from to 
those from ^4° exceeds 3 to 1. The dominance in this zone stems from the H'^-x'^-xl 
coupling {H^-XtXs coupling) being two to three times larger (smaller) than the A'^-Xs-Xa 
coupling (A°-X3-X3 coupling), combined with the fact that the xtxl decays are about twice 
as likely to produce Ai events as those of X^Xs- This of course means that Xa, has a higher 
leptonic BR than Xs- This in turn is due to X3 decaying into Xi^^ about half the time {Z^ 
gives lepton pairs ~7% of the time), while xt almost never decays this way, instead having 
larger BRs to charged sleptons [and xf VF''^] which always [~21% of the time] yield charged 
lepton pairs. The situation changes quickly once the H^, A^ XaXa^ X2X2 thresholds are 
(almost simultaneously, see Fig. 8) crossed, thereafter for higher Ma values the ^4" and 
contributions remain reasonably close to each other as already stated. 

As with points in Fig. 7, direct chargino/neutralino pair- production and slepton pair- 
production together with Z^^*'> Z^^*'> production make up most of the background surviving 
the cuts. Now, however, these are joined by a minor segment due to tH~ -\- c.c. production, 
which depends on Ma in addition to tan (3. 

Results showed gh — > tH~ + c.c. could yield several events at points in the discovery 
region. Since the presence of a charged Higgs boson would also signal that there is an extended 
Higgs sector, these events could easily have been grouped with the signal rather than with 
the backgrounds. Clearly though the set of cuts used in this work is not designed to pick 
out such events. The jet cut typically removes roughly two-thirds to three-quarters of these 

^'''Hcrc are some results from specific points in this region: for Ma = 510 GcV and tan /3 = 10, 16, 25, 40, the 
percentage of A° signal events (again, after cuts, excluding the M inv. m. cut), is 23%, 17%, 11.5%, 21%. 
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events. Here though it is interesting to note that, despite the presence of a top quark, the 
jet cut does not remove all such events (unlike results found for squark and gluino events 
and four-lepton ttX events). A more effective set of cuts for tH~ , tH'^ events is developed 
in [20], wherein substantially larger numbers of charged Higgs boson events survive the cuts 
therein at favorable points in the MSSM parameter space. It is also worth noting though 
that the reach of the discovery region (at a favorable point in the MSSM parameter space) for 
the H^,A^ — > 4£ signal as described in this work surpasses that of the charged Higgs boson 
discovery regions found in [20]. (or in any other previous work on Higgs boson decays to 
sparticles) . 

An aspect to be mentioned in this connection, already highlighted in Ref. [24], is the 
somewhat poor efficiency for the signals following the Z°-veto, especially when combined with 
the fact that the background survives the same constraint. On the one hand, a 

non-neghgible number of events in the signal decay chains leading to 4^A'" final states actually 
proceed via (nearly) on- mass-shell Z° bosons, particularly for MSSM Point 2, in which the 
mass differences m~o — m~o {i = 3, 4) can be very large, unlike the case m~o — m-o for MSSM 
Point 1 (and in previous studies limited to only X2X2 decay modes). On the other hand, 
the rather large intrinsic Z° width (when compared to the experimental resolution expected 
for di-lepton invariant masses) combined with a substantial production cross-section implies 
that Z°^*^Z^^*^ events will not be totally rejected by the Z°-veto. Altogether, though, the 
suppression is much more dramatic for the Z^Z^ background than for the signal, and so this 
cut is retained (though the Z^-veto will be dropped in some instances in the context of the 
forthcoming wedgebox analysis). Also, varying the size of the 10 GeV window around Mz did 
not improve the effectiveness of this cut. 

4.2 mSUGRA benchmark points 

Turning attention briefiy to the results within the more restrictive mSUGRA framework for 
SUSY-breaking, results for mSUGRA Point A and mSUGRA Point B (as defined in Sect. 3) 

are presented in Tables 5-6. Mass spectra for these parameter sets arc given in Table 1. For 
mSUGRA Point A ample signal events are produced and survive the cuts to claim observation 
of the Higgs boson at 100 fb~^. The largest background is from direct slepton production, with 
direct neutralino/chargino production also contributing significantly, whereas SM backgrounds 
are virtually nil. Note how the cap suffices to efiminate the background from colored 
sparticle (squarks and gluinos) production. 

Recall that for mSUGRA Point A the signal is dominated by H^, X2X2 decays, 

whereas for mSUGRA Point B heavier -inos make major contributions. Thus, a wedgebox 
plot analysis of the former is expected to show a simple box topology, while in the case of the 
latter there unfortunately may be too few events (even with 300 fb~^ of integrated luminosity) 
to clearly discern a pattern. For mSUGRA Point B, 9(10) signal events survive after all cuts 
(save the 4£ inv. m. cut), while 6 background events survive, assuming 100 fb^^ of integrated 
luminosity. This is insufficient to claim a discovery by the criterion of Relation (2). However, 
when the integrated luminosity is increased to 300 fb^^, then the raw number of signal events 
suffices to cross the discovery threshold. Unfortunately though, for mSUGRA Point B the 
background from colored sparticle production is not removed by the upper limit imposed on 
One can however stiffen the E'^?* cut, capping the allowable jet transverse energy at 
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30 GeV rather than 50 GeV and thus ehminate much of this background without diminishing 
the signal rate significantly. Then, with 300 fb~^ of integrated luminosity the discovery criteria 
can be met. 

An earlier ATLAS study [47, 16] also sought to map out the discovery reach of the Higgs 
boson to neutrahno four-lepton signature within the mSUGRA framework transposed onto the 
{Ma, tan /3) plane. Though some statements to the contrary are included in this ATLAS study, 
it does seem to have been focused on the X2X2 contributions (analogous to previously-discussed 
general MSSM studies of this signature), thus apparently omitting parameter sets such as 
mSUGRA Point B considered herein. Thus, the viability of mSUGRA Point B indicates 
an enlargement of the signal discovery region to higher values of Ma (and the mSUGRA 
parameter Mq) at intermediate values of tan/5 (i.e., in the 'decoupling' region) from that 
reported in this ATLAS study (akin to the enlargements shown in the general MSSM case, 
though the extent of this enlargement in the case of mSUGRA models will not be quantified 
herein) . 



Table 5: Event rates after the successive cuts defined in the text for mSUGRA Point A 
(assuming an integrated luminosity of 100 fb~^). 
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Table 6: Event rates after the successive cuts defined in the text for mSUGRA Point B 
(assuming an integrated luminosity of 100 fb~^). 
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5 Wedgebox analysis of Higgs boson decays to — ino 
pairs 

The wedgebox plot technique was introduced in a previous work [27] which focused on neu- 
trahno pairs produced via colored sparticle production and subsequent 'cascade' decays. An- 
other work [48] has just recently focused on neutralino pairs produced via EW processes, 
includiner via a boson or via H^,A^ production; the former is termed 'direct' production 
while the latter is 'Higgs-mediated' production. A jet cut was found to be fairly efficient in 
separating these two neutralino pair-production modes from cascade production assuming the 
colored gluinos and squarks are fairly heavy. 

To utilize the wedgebox technique, the criteria for the final four-lepton state arc further 
sharpened by demanding that the final state consist of one e+e" pair and one fi^fi^ pair 
The wedgebox plot then consists of the M{fj,^fj,^) invariant mass plotted versus the M(e+e^) 
invariant mass for all candidate events. If a given neutralino, x°, decays to the LSP, x?, and a 
charged lepton pair via a three-body decay mediated by a virtual or virtual/off-mass-shell 
charged slepton, then M(£'^£~) is bounded from above by m~o — m~o (and from below by 
if lepton masses are neglected). Given a sufficient number of events, the wedgebox plot 
of the signal events will be composed of a superposition of 'boxes' and 'wedges' [27], in the 
M(e'^e~)-M{ij,~^ H") plane resulting from decay chains of the form: 

i^°,^°-X°X;-e+e>Vx?X?. (3) 

If x° decays into an e~^e~ {^^ ji") pair, then M(e+e~) (M(/x^yU~)) is bounded above by 
m~o — m-o (m~o — ""^^o)- On the other hand, if (x^) decays into a ji'^ ji~ (e+e^) pair, 
then these M(e^e~) and M{fi~^fi~) upper bounds are swapped. Superposition of these two 
possibilities yields a 'box' when i = j (which will be called an H-i box') and a 'wedge' (or 
'L-shape') when i j (this will be called an 'i-j-wedge'). 

A heavy neutralino, ; could instead decay to the Xi + leptons final state via a pair of 
two-body decays featuring an on-mass-shell charged slepton of mass^^ mj. Events containing 
such decays will lead to the same wedgebox pattern topologies as noted above; however, the 
upper bound on M{i~^i~) is modified to 




M{i+i-) < m~o, 1 - I z:^ K 1 - I TT^ J ■ (4) 



The M(£+£~) spectrum is basically triangular in this case and sharply peaked toward the upper 
bound, while the former three-body decays yield a similar but less sharply peaked spectrum. 
The two-body decay series alternatively could be via an on- mass-shell resulting in an 
M{f+£-) = Mz spike. 

Additional complications can arise if the heavy neutralino Xi can decay into another neu- 
tralino x'i (j 7^ 1) or a chargino which subsequently decays to yield the Xi ^T^sd state. These 



^^In fact, this extra restriction is not strictly necessary, since recent preliminary work shows same-flavor 
four-lepton final states can be correctly paired with a reasonably high efficiency for at least some processes 
and some points in the MSSM parameter space [49] . 

^''Note that this is the physical slepton mass, not the soft mass input. 
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may introduce new features to the wedgebox plot: Xi to Xj (j 1) decay chains involving 
xl — > i^i~X2: Xa ~^ ^'^^~X2^ and/or x\ ^"^^"Xs will generate additional abrupt event 
population changes or edges, termed 'stripes,' on the wedgebox plot. One can imagine quite 
elaborate decay chains, with X4 — > X3 ^ X2 ^ Xi fo^ instance. However, such elaborate chains 
are very unlikely to emerge from any reasonable or even allowed choice of MSSM input param- 
eters. Further, each step in such elaborate decay chains either produces extra visible particles 
in the final state or one must pay the price of the BR to neutrino-containing states. The latter 
tends to make the contribution from such channels insignificant, while the former, in addition 
to also being suppressed by the additional BRs, may also be cut if extra restrictions are placed 
on the final state composition in addition to demanding an e+e" pair and a yU^/i" pair. The 
aforementioned extra visible particles could be two more leptons, meaning that all four leptons 
come from only one of the initial -inos, x° — > (■'^(■~Xk ~^ i'^i''i'~^i'~Xi: while the other -ino, 
which must yield no leptons (or other visible final state SM particles forbidden by additional 
cuts), decays via x'j ~^ ^^x\ or QQXv Again though such channels will be suppressed 

by the additional required BRs. A further caveat is that decays with extra missing energy 
(carried off by neutrinos, for example) or missed particles can further smear the endpoint. 
The presence of charginos may also further complicate the wedgebox picture. Heavier -inos 
can decay to the LSP -|- lepton pair final state via a chargino, Xi ~^ ^^1^X1 ~^ i'^ui''~P'xi: 
or a Higgs boson itself may decay into a chargino pair, with one chargino subsequently yield- 
ing three leptons while the other chargino yields the remaining one (such events are called 
'3-1-1 events' [48]). The chargino yielding three leptons will typically decay via a X2^ resulting 
in a re-enforcement of the solely -ino-generated wedgebox plot topology. A single chargino- 
generated lepton paired with another lepton from a different source produces a wedge-like 
structure but with no definite upper bound. For a more in-depth discussion of these nuances, 
see [48]. 

The right-hand plot in Fig. 9 shows the wedgebox plot obtained in the case of MSSM Point 
1, assuming an integrated LHC luminosity of 300 fb~^. Criteria for event selection are as given 
in the previous section, save that the more restrictive demand of an e'^e~ ^'^ix~ final state is 
applied while the Z^-Yeio and four-lcpton invariant mass cuts are not applied. Both signal 
and background events are included; the former are colored black. The latter consist of both 
SM backgrounds (on- or off-mass-shell Z°-boson pair-production — z'^^*^Z^^*\ 83 events, and 
ttZ'^^*\ largely removed by the missing energy and jet cuts, 2 remaining events; these events 
are colored red and purple, respectively, in Fig. 9 ) and MSSM sparticle production processes 
('direct' neutralino or chargino production, 4 events, and slepton pair-production, 22 events; 
such events are colored green and blue, respectively, in Fig. 9). No events from colored sparticle 
production survive the cuts, particularly the jet cut — this is a crucial result. Signal events 
consist of 14 events and 25 events, yielding a signal to background of 39 : 111 = 1 : 2.85. 
With S/y/B — 3.7, this is not good enough to claim a discovery based on Relation (2). If the 
input CP-odd Higgs boson mass is lowered to Ma = 400 GeV, whose wedgebox plot is the left- 
hand plot of Fig. 9, then the number of signal events rises to 14 + 52 = 66 and A° events 
(runs for MSSM backgrounds gave 2 'direct' neutralino-chargino events and 26 slepton-pair 
production events), yielding S/\/B = 6.2 and satisfying Relation (2). Note how the increase 
is solely due to more 74°-generated events. Comparing the Ma — 500 GeV (MSSM Point 1) 
plot (b) and the Ma = 400 GeV plot (a) in Fig. 9 shows how the increased number of signal 
events in (a) more fully fills in the 2-2 box whose outer edges (dashed lines in the figure) are 
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Figure 9: Wedgebox plot for MSSM Point 1 inputs with Ma = 500 GeV (6) and shifting to 
M4 = 400 GeV (a), assuming an integrated luminosity of 300 fb~^. Neither the Z'^-veto cut 
nor the 4-lepton invariant mass cut are enabled. Black-colored markers are for the and 
A'^ signal events. SM background events from Z'^^*^Z^^*^ (where either one or both of the Z°'s 
are permitted to be off-mass- shell are red), while the two surviving tiZ^^*^ events are purple. 
MSSM background events from slepton production or direct neutralino/chargino production 
are in blue and green, respectively. The horizontal and vertical dashed lines forming a box 
are at the location m~o — m~o. MSSM Point 1 -ino and slepton inputs are fi = —500 GeV, 
M2 = 180 GeV, Ml = 90 GeV,' mj = m~ = 250 GeV. 

given by m~o — m~o = 86.6 GeV since for these input parameters slepton masses are too high 
to permit X2 decays into on-mass-shell sleptons. 

A key observation is that the distributions of the signal and the background events differ 
markedly^^. All but one of the signal events lie within the 2-2 box^^. The majority of the 
slepton pair-production events (19 out of 26 events for (a) and 17 out of 22 events for (6)), 
the dominant MSSM background, lie outside the 2-2 box. The topology of these '3+1' events 
is a 2-2 box plus a wedge lacking a clear outer edge extending from said box (see [48]). The 
few 'direct' neutralino and chargino production events happen to all lie within the 2-2 box; 
however, these events are actually due to^° X2X3 pair-production and thus, for a larger sample, 

^®0n the other hand, the distributions of A'-* and events show no substantial systematic differences in 
their distributions' wedgebox plot topologies. 

^^Note that a similar result is found in Fig. 16 of [9]. There, however, only signal events were shown, 
and, since a priori only H'^^A'^ —f X2X2 events were considered, the vast array of other potential wedgebox 
topologies was not brought to light. 

■^^If direct neutralino pair-production produces a significant number of events, then the dominant source of 
said events is always X2X3 production; X2X2 production is heavily suppressed. See discussion in [48]. This 
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such events would populate a 2-3 wedge with many of the events falling outside of the 2-2 
box. 

SM background events are concentrated on and around lines where either M(e+e~) and/or 
M[ix^ ix~) equals M^, which unfortunately is close to the outer edges of the 2-2 box. Using 
the unfair advantage of color-coded events, one can correctly choose to place the edges of 
the box so as to exclude most of the SM background events. Experimentalists may have a 
more difficult time deciding on wedgebox edges that lie too close to Mz- Though, at the 
price of perhaps losing some of the signal events^^, one could make a selection rule of an 
effective 2-2 box with edges sufficiently within Mz in such cases. Correct identification of 
the outer edge value for the 2-2 box removes all but 11 of the 85 SM background events. 
The signahbackground is then 39 : 20 for (6) and 66 : 19 for (a), an immense improvement 
in the purity of the samples — both points now certainly satisfy the Relation (2) criterion. 
Accepting only points lying within a box with outer edges at 80 GeV, more safely eliminating 
SM events, leads to a signahbackground of 33 : 12 for (6) and 59 : 14 for (a). Note 

that one can also select points lying well outside the 2-2 box to get a fairly pure sample (at 
this point in the parameter space) of slepton pair-production events. Even if one does not 
know where Nature has chosen to reside in the MSSM input parameter space, the selection 
of only events occupying one distinct topological feature of the experimental wedgebox plot 
may yield a sample pure enough (though one may not know exactly what purified sample one 
has obtained!) to be amenable to other means of analysis (perhaps entailing some addition 
reasonable hypotheses as to what sparticles might be involved) [51]. 

Fig. 10 in turn examines several related choices for input parameter sets, including MSSM 
Point 2 — which is plot (c) therein, in which and have large BRs into heavier -ino pairs 
such that the majority of the M signal events do not arise from X2X2 decays for all points save 
that of plot (6). Plot {d) differs from MSSM Point 2, plot (c), only in that the Higgsino mixing 
parameter fi is shifted from ji = —200 GeV to /i = —250 GeV. Yet even this modest change 
drastically alters the topology of the resulting wedgebox plot. This is illustrative of how the 
wedgebox plot may be useful in extracting fairly detailed information about the -ino spectrum 
and corresponding MSSM input parameters. In plots (a) and (b) of Fig. 10 the EW gaugino 
input parameters are raised from M2 = 200 GeV in plots (c) and {d) to M2 = 280 GeV (recall 
the assumption used herein that the value of Mi is tied to that of M2). Also tan/3 is lowered 
from 35 to 20, while // values of plots (c) and (d) are retained. Again, these shifts in input 
parameters radically alter the resulting wedgebox topology. Plots (a) and (6) clearly show 
wcdgc-likc topologies. Note again the markedly different event distributions for the signal and 
background events in all four plots, but particularly striking in plot (a). Note how the four 
MSSM parameter set points yielding the wedgebox plots depicted in Fig. 10 all might crudely 
be categorized as high tan/3, low low to moderate M2, and light slepton points. Yet the 
associated wedgebox plots come out decidedly different. 

Taking advantage of knowing which points in MSSM parameter space are being simulated 

leads to the general conclusion that, with a jet cut in place to remove cascade events from colored sparticle 

decays, the appearance of a disproportionately strong (densely populated) box on a wedgebox plot is highly 
indicative of the presence of Higgs-boson-generated events. The caveat to this being that chargino production 
can generate a box-shape in some rather limited regions of the MSSM input parameter space. Again, see [48] 
for further discussion. 

Correct edge values from which to reconstruct information on the -ino mass spectrum would also be lost. 
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Figure 10: Wedgebox plot for MSSM Point 2 inputs (c) and shifting to Ma = 400 GeV 
(left), assuming an integrated luminosity of 300 fb^^. Neither the Z°-veto cut nor the 4- 
lepton invariant mass cut are enabled. Black-colored markers are for the and signal 
events. SM background events from Z^^*^Z^^*^ (where either one or both of the Z°'s are 
permitted to be off-mass- shell) are red, while the two surviving tiZ^^*^ events are purple. 
MSSM background events from slepton production or direct neutralino/chargino production 
are in blue and green, respectively. The horizontal and vertical dashed lines forming a box 
are at the location M~o — M~o. MSSM Point 1 -ino and slepton inputs are /i = —500 GeV, 



M2 = 180 GeV, Ml = 90 GeV, mj = m~ = 250 GeV. Also indicated by dashed lines on 

^soft soft 

the plot are the 2-2, 3-3 and 4-4 box edges found from relation (4) 
edges for (c) which are from m~o — '"^5^0 • 



save for the 2-2 box 
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Table 7: Percentage contributions to H^,A^ 4i events from the various neutralino and chargino 
pair-production modes for the four MSSM Parameter set points given in Fig. 10. Based upon 
IS A JET (IS ASUS Y) 7.58 [30] with no consideration given to any cuts. Decays that are kinematically 
not allowed are marked by a 0; contributions below 0.001% are marked as negligible (neg). H^, 
^o(*)^o(*)^ jjQ _^ ^0^0 g^^^ ^0 ^ /i0^o(*) j^aj^g negligible contributions in ah cases. Also given are 

the number of H^jA'^ signal events and the number of background events, assuming 300 fb-^ of 
integrated luminosity as in the figure. 
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(something the experimentalist cannot know in the actual experiment) allows comparison 
between the assorted calculated production rates at the four points and the observed features 
on the wedgebox plots. Table 7 gives such theoretical estimates based on analysis of ISAJET 
(ISASUSY) 7.58 results for the four points^^. It must be borne in mind though that effects 
from cuts may alter the percentage contributions found on the wedgebox plots from those 
given in Table 7. 

The first thing to notice from this table is the virtual absence of events stemming from xl 
to x5 decays for MSSM Point 2 = plot (c) relative to the other three points. This is due to 
the fact that, for this input parameter set, the sparticle spectrum satisfies the condition that 
m~ < m~o < mj^, meaning that xl mainly decays via an on-mass-shell sneutrino 'spoiler' 
mode, xl ^ ^ Xi^^j si'iid its BR into a pair of charged Icptons is highly suppressed. For 
the other three points, m~o > mjj.,m~. Actually, of the four wedgebox plots shown in Fig. 10, 
the one for MSSM Point 2 most closely resembles a simple box. However, Table 7 indicates 
that (before cuts) 45.0% of the events are from X3X4, 39.6% of the events are from xlx% 
12.7% of the events are from X1X21 X2X2 ■ 

In Fig. 10, charged slcptons arc now light enough so that the neutralino to slepton decay 
chains, which make significant contributions to the four-lepton signal events, may proceed 
via on-mass-shell charged sleptons. So while the outer edges of the 2-2 box in Fig. 9 was 

^^Table 3 given previously corresponds to column (c) in Table 7 with the and contributions listed 
separately. 
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determined by the X2~Xi mass difference, here Relation 4 brings the slepton masses into 
play^^. In plot (a), virtually all Xi to Xi decays proceed via on- mass-shell sleptons, but 
only the xi to Xi decay edge is significantly altered (by more than a couple GeV) — from 
m~o — m-o = 185 GeV to 151-156 GeV (at this point, 18% of four-lepton events are from XsXa 
according to Table 7). On the other hand, in plot (6), where the Xa ^-Iso decay to Xi via 
on-mass-shell sleptons, edges are shifted from m~o — m-o = 82, 124, 192 GeV to 76-78, 101- 
107, 140-149 GeV for i = 2,3, 4, respectively^'', with i = 2, 3, 4 decays all making noteworthy 
four-lepton event contributions. For MSSM Point 2 = plot (c) , the shift in the X3 to Xi decay 
edge is only 3.5-5 GeV while the Xa to Xi edge is virtually unchanged. This accounts for 87.3% 
of the four-lepton events by Table 7. The situation with X2 slightly complicated: X2 '^^^ 
only decay into Xi via an on-mass-shelP^ Jli, and this would lead to a tremendous shift in the 
edge position (from 61 GeV to 15 GeV); however, this is so close to the kinematical limit that 
decays through off-mass-shell should be competitive (again placing the edge at ~61 GeV). 
But, since X2 decays lead to only a tiny fraction of the four-lepton events, note how there 
is no visible edge or population discontinuity at this location (the innermost dashed box) on 
the wedgebox plot. Lastly, with plot (d), again on-mass-shell slepton decays totally dominate 
for i = 2,3,4, but only the X2 to X? decay edge is significantly shifted (from 75.2 GeV to 
51.5-60.4 GeV^^. But, by Table 7, this decay is the most important contributor to the signal 
events. 

For plot (a) of Fig. 10, the expected 2-4 wedge stands out clearly among the signal events, 
with outer edges at the expected location. The background is mostly from direct X2X3 direct 
production, giving the 2-3 wedge shown in green (direct neutralino-neutralino production is 
predominantly X2X3 ^t all interesting points in the MSSM parameter space, with direct X2X2 
production always highly suppressed [48]). The proximity of this wedge's outer edges to the red 
Mz lines may complicate the experimental analysis; however, if the SM Z^^*^Z^^*^ background 
is well-modeled, a subtraction technique to clear up this zone may be feasible. Note that 
selecting only events with 100 GeV < M(e+e-) < 150 GeV, < M(/x+/x-) < 50 GeV or 
< M{e+e ) < 50 GeV 100 GeV < M{ii+ii-) < 150 GeV, corresponding to the legs of the 
2-4 wedge lying beyond the 2-3 wedge and the Z°-line, changes the signalrbackground ratio 
from 728:683 seen on the plot to 128:15. This is an example of a cut that can be applied a 
posteriori based on the examination of the wedgebox plot — as opposed to assuming a priori 
extra knowledge about where in the MSSM parameter space Nature has chosen to sit. 

Plot (b) of Fig. 10 mainly shows a densely-populated 2-2 box whose edges are well inside 
the Mz hues. A faint 2-3 or 2-4 wedge is also discernible (in fact Table 7 shows this to be a 

^•^Unfortunately, the physical slepton masses input into HERWIG 6.5 are generated by ISASUSY 7.58 [30], 
which neglects a left-right mixing term cx m|/x^ tan^ /3 (see [27]). While this term is negligible for selections, it 
does shift the physical smuon masses by as much as a few GeV. Neglecting this term results in degenerate soft 
slepton inputs leading to degenerate physical selectron and smuons masses (so the smiion masses for MSSM 
Point 2 given in Table 1 are changed into the mass values given there for the selectrons), which in turn may 
noticeably under-estimate the mass splitting between smuons and thus the thickness of the edges shown on 
the plots. Later versions of ISAJET correct this oversight, as do private codes employed in Sect. 2. 

■'"'Due to the program oversight mentioned in the last footnote, the thicknesses of these edges shrink to 
75.7-76.5, 103.4-104.9, 143.5-145.9 GeV, respectively. These values are represented by the dotted lines on the 
plots. 

3^ Again, this feature is lost in HERWIG 6.5/ISAJET 7.58 . 

36ln HERWIG 6.5/ISAJET 7.58 this width shrinks to 55.6-57.0 GeV. 
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2- 4 wedge), while the empty upper-right corner which does not join with the 2-2 box suggests 
that X2X4 ^^'^ XsXa decays are present while xlxl absent (further suggesting that said 
decay mode is kinematically inaccessible, which helps pin down the relative masses of the 
heavy Higgs bosons and the heavier neutralinos). 

Plot (c)'s most obvious feature is an outer box, which in fact is a 4-4 box. Topology alone 
does not distinguish this from a plot dominated by a 3-3 box or a 2-2 box, though the location 
of the outer edges well beyond Mz might give pause for entertaining the latter possibility. A 

3- 4 wedge may also be discerned from the somewhat diminished event population in the upper 
right-hand box in the plot. Comparison of this plot with the other three quickly points out the 
absence of a dense event-population in this plot. Seeing such a wedgcbox plot experimentally 
strongly hints that leptonic X2 decays are being suppressed, perhaps with a mass spectrum 
favoring sneutrino spoiler modes as noted above. 

Like plot (6), plot (d) shows a 2-2 box, but with outer edges at a very different location. 
Plot (d) also has more signal events outside of the 2-2 box than does plot (b), and said events 
are more scattered in (d). A lot of these events are from /J", A" decays into xfx2 pairs. Thus, 
the alignment of the wedgebox features to the dashed lines derived from neutralino features 
shown is less compelling. 

In both Fig. 9 and Fig. 10, note how closely the wedgebox plot features, obtained by the full 
event generator & detector simulation analysis, conform to the dashed-line borders expected 
from the simple formula 4. This strongly supports the assertion that a wedgebox-style analysis 
is realistic in the actual experimental situation. 

6 Summary and conclusions 

Recapping the findings presented herein: 
6.1 New signals 

For many interesting choices of the basic input parameters of the MSSM, heavier Higgs boson 
decay modes of the type H^,A^ — > XiX^: with i,j 7^ 1 are potentially important LHC signal 
modes. The neutralinos' subsequent leptonic decays, typified by x^ can yield a 

four-isolated-lepton (where here i refers to electrons and/or muons) plus missing-transverse- 
energy signature. Such leptonic neutralino decays may proceed via either an intermediate 
charged slepton or via an intermediate Z^^*\ where in either case this intermediate state 
may be on- or off-mass-shell. The present study presents for the first time a systematic 
investigation of the potential for discovering such a signature at the LHC, including all possible 
such neutralino pairs: XsX". X2X3> X2X4> andx4X4- Other Higgs boson decays that 

may lead to the same signature are also incorporated, inchiding: decays to chargino pairs 
H^,A^ — > XiXS 1X2X2 1 iri which case xS yields three leptons while the other chargino gives 
the fourth; H^,A^ XiX3,XiX4) where the X3 or X4 must provide all four leptons; and 
//o ^ AO ^ hPZ^^*\ & H^,A^ I+I-, all three of which yield negligible 

contributions in all cases studied. This surpasses previous studies which restricted virtually 
all of their attention to H^, A^ X2X21 ^^^o did not consider the possibility of neutralino 
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decays to on-mass-shell sleptons (with the incorporation of the heaviest neutralinos as is done 
herein this assumption becomes particularly restrictive). 

Naturally, at least some of the -inos must be reasonably hght for this H^,A^ — > 4:i + E^^^^ 
signature to be seen. Parameter-space scans studying the potential scope of such a signal 
indicate that the -ino parameter M2 needs to be relatively low while the Higgsino mixing 
parameter /j, need not be so constrained (however, if is not also relatively low, then the 
signal is dominated by the X2X2 mode). Relatively light slepton masses are also quite helpful, 
and the slepton mass spectrum plays a crucial role in determining for what values of the other 
MSSM input parameters large rates may occur. Said large rates are possible throughout 
most of the phenomenologically-intcrcsting value ranges of the Higgs-sector parameters Ma 
and tan /?, depending of course on the accompanying choice of other MSSM inputs, as the 
discovery regions delineated herein illustrate. 

6.2 Comparison with previous results 

To clearly demonstrate the potential importance of the H^, — > 4£ -|- E^^^^ signature in the 
hunt for the heavier Higgs bosons. Figs. 11 and 12 again show the discovery regions associated 
with MSSM Point 1 and MSSM Point 2 neutralino input parameter sets (as depicted before in 
Figs. 7 and 8, respectively), but this time with a logarithmic scale for tan/? and also showing 
the expected reaches, assuming 300 fb~^ of integrated luminosity at the LHC, of Higgs boson 
decay modes into SM daughter particles as developed by the ATLAS collaboration [45]^^. 
Clearly, the new neutralino decay mode signature can extend the discovery reach for the 
heavier MSSM Higgs bosons to much higher values of Ma, and also offer at least partial 
coverage of the so-called 'decoupling region' where only the lightest Higgs state could be 
established in the past (through its decays into SM objects) and where said may be difficult 
to distinguish from the sole Higgs boson of the minimal SM. Thus, a more complete analysis 
of the A° — > x°x° modes as is presented here may be crucial to the establishment of an 
extended Higgs sector. The inclusion of the heavier neutralinos, Xs xl^ absent in previous 
studies, is essential in extending the reach of the H^jA'^ A£ + E^^^^ signature up to the 
higher Higgs boson masses unattainable by the SM decay modes. 

It should be noted that the ATLAS discovery contours presented in Figs. 11 and 12 are not 
obtained using the same choice of MSSM input parameters as are the H^,A^ — > x°x° discovery 
regions developed in the present work. In fact, the ATLAS discovery regions used input choices 
designed to eliminate, or at least minimize, the Higgs boson decays into sparticles. Thus, the 
reach of the ATLAS discovery contours essentially represents the maximum expanse in the 
MSSM parameter space achievable through these Higgs boson decays to SM particles under 
the (unsubstantiated) assumption of a very heavy sparticle sector. Stated another way: were 
the ATLAS discovery regions to be generated for the same set of neutralino input parameters 
as the H^, A^ XaXj discovery regions presented herein, the former may well shrink in size 
(and certainly not increase), further emphasizing the importance of thoroughly studying the 
H^,A^ 4:£ + E^'^^^ signature. It would certainly be desirable to re-do the SM-like signature 
reaches of MSSM Higgs bosons in the presence of light sparticle spectra identical to those 



ATLAS collaboration discovery region contour lines in Figs. 11 and 12 have been remade to match as 
closely as possible those in the original plot. 
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Figure 11: Discovery regions in the (M^jtan/?) plane, here with a logarithmic tan/3 scale, 
assuming MSSM Parameter Set 1 -ino inputs and for Cint = 100 fb~^ and 300 fb~^, for the 
(lower plot) MSSM Higgs bosons' 4i signals from their decays into neutralino or chargino pairs 
(here H^,A^ decays to X2X2 totally dominate). This is shown juxtaposed (upper plot) with 
300 fb~^ regions for MSSM Higgs boson signatures from decays to SM particles based upon 
LEP results and ATLAS simulations [45], where labels represent: 1. Z^Z^* 4 leptons; 

2. t_-^ bH+, H+ T+u + c.c; 3. bb; 4. h° 77 and W^hytth^, h° -> 77; 

5. bbH^,bbA^ with H^A^ bb; 6. H+ ^ tb + c.c; 7. H^[A^ fi+fi-; 8. H^A^ t+t-_; 
9. gb iH+, H+ r+u + c.c; 10. H° h^h^ bhr, H- Z^h^ -> i^i~bb; 

12. H^/A^ -> tt. Note that SM discovery regions are not for the same input parameters: 
they presume a very heavy sparticle spectrum; identical MSSM inputs to those used for the 
lower plot may well yield smaller SM discovery regions in a revised upper plot. For the 4i 
signals from x'iXpXtiXn decays, the MSSM Parameter Set 1 -ino/slepton parameters are 
u = -500 GeV, M2 = 180 GeV, Mi = 90 GeV and mj = m~ = 250 GeV. 

'^soft 'soft 
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Figure 12: Discovery regions in the [Ma, tan (3) plane, here with a logarithmic tan (3 scale, 
assuming MSSM Parameter Set 2 -ino inputs and for Cint = 100 fb~^ and 300 fb~^, for the 
(lower plot) MSSM Higgs bosons' Ai signals from their decays into neutralino or chargino 
pairs (here Higgs boson decays to higher-mass neutralinos typically dominate). This is shown 
juxtaposed (upper plot) with 300 fb~^ regions for MSSM Higgs boson signatures from decays 
to SM particles as in Fig. 11. For the M signals from XiXp xtiXn decays, the MSSM Parameter 
Set 2 -ino/slepton parameters are u = -200 GeV, = 200 GeV, Mi = 100 GeV, mj = 

^soft 

150 GeV and m~ ^ = 250 GeV. Here Higgs boson decays to a variety of higher mass -inos 
(see text) constitute the majority of the signal events. Note that, as in Fig. 11, since ATLAS 
discovery regions presume a very heavy sparticle spectrum, SM discovery regions made for the 
same MSSM input parameters as used in the lower plot may well yield smaller SM discovery 
regions in a revised upper plot. 
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studied herein for the Higgs-to-sparticle decay channels; however, this is clearly beyond the 
scope and capabilities of this study. It also must be emphasized that the diminution of the 
expected signatures from SM decay modes of the MSSM Higgs bosons was investigated in [7] 
and thus is fairly well-established as well as inherently sensible. 

Previous studies exploring Higgs-to-sparticle decay channels, whether for neutral Higgs 
bosons {e.g., CMS [9]) or for charged Higgs bosons {e.g., ATLAS [52], CMS [20]), — and 
comparing, to some extent, SM and SUSY decay modes — have not re-scaled the reaches 
of previously-studied SM decay channels (done by the same collaboration) to allow a rea- 
sonable comparison to the new-found sparticle decay modes; nor have the SM decay modes 
been re-analyzed for the same set of MSSM input parameters. Yet clearly such comparisons 
are absolutely essential to gauge the scope and impact of the new sparticle-decay channels. 
Certainly, the comparisons presented in Figs. 11 and 12 are less than optimal; however, they 
are far from un-informative. 

It is also important to keep in mind that the assumptions inherent in the ATLAS (and 
CMS) discovery regions for the SM decay modes of the MSSM Higgs bosons are no less 
restrictive than the choices of MSSM input parameters made to generate the two M + E"™®^ 
discovery regions in this study. The parameter space scans of Sect. 2 further enable the reader 
to put the two discovery regions shown here into a wider perspective. 

6.3 Production and decay phenomenology of the signal 

The new H°,A° — > + E^'^ discovery regions have been mapped out using a full event 
generator-level analysis utilizing HERWIG coupled with a detector simulation on a par with 
experimental analyses. All significant backgrounds have been included in the analysis, some 
for the first time in the study of such a signature. The importance of the restriction on 
jet activity employed herein is particularly noteworthy. Without such a cut the Higgs signal 
could be swamped by the cascade decays of colored sparticles (gluinos and squarks) , unless said 
sparticles are a priori assumed to be quite heavy (at or above the TeV scale). The ultimate 
limit of this type of jet cut, to demand that events be 'hadronically quiet' quickly springs 
to mind as an attractive search category. Yet care must be taken here since, in Higgs boson 
production via gg — > A° and bb — > jets emerge in the final state alongside the Higgs 

bosons due to PS effects, though such additional jets tend to be rather soft and coUinear to the 
beam directions. In addition, rather than emulating Higgs boson production via gg — > ^4° 
and bb H^, one could instead consider gg ggH^, ggA^ and gg — > bbH^, U}A^ processes, 
in which case one might worry about stronger jet activity emerging. The true signal rate is 
the sum of these and the previous process types, after making a correction for the overlap 
(as discussed previously). HERWIG simulations of gg — > bbH^,bbAP at selected points in 
the parameter space indicate that the these processes are in fact removed by the jet cut 
imposed herein. To better optimize the level of hadronic activity that should be allowed, full 
implementation of 2 3 loop processes {gg — >■ ggH^, ggA^ and other channels yielding two 
light jets and a H^, A^ in the final state) into HERWIG must be completed (work in progress 
[53]). 

The BRs of and ^4° to the assorted -ino pairs can certainly differ markedly in regions 
where the signal is large, as seen for instance in Tabic 3; thus one must not assume that 
the two contribute a roughly equal number of events to the 4i + E^^^^ signal rate. On the 
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other hand, results also show that only in quite narrow low-M^i threshold regions within the 
discovery areas (wherein the small M^-Ma mass difference is crucial) do events due to one or 
the other Higgs boson (in this case the lighter ^4°) totally dominate, producing in excess of 
90% of the signal events. General statements beyond this concerning the and A° admixture 
present in the signal seem elusive. Throughout the X2X2"dominated discovery region of Fig. 
7, produced the majority of the events (though in some cases only slightly more than 
whereas in Fig. 8 there were substantial zones in which events dominated (as well as 
large segments wherein the two Higgs boson contributions were within ~20% of each other). 
Finally, though the cuts did typically eliminate slightly more events than events, this 
effect was of little significance. 

6.4 The topology of the signals 

Note that in comparing the signal with the MSSM backgrounds, the present study follows 
the standard procedure of comparing signal and background rates at the same point in the 
MSSM parameter space. One could well ask whether or not larger backgrounds at a different 
point in parameter space could lead to the number of excess events attributed to the signal 
at the designated point in the MSSM parameter space. One way of addressing this issue 
is to look at the distribution of the signal+background events on a M(e"'"e~) vs. M{fi'^^~) 
wcdgcbox plot in addition to merely asking what is the raw rate. To wit, analyses of selected 
points in parameter space, again at the full event generator + detector simulation level, are 
presented illustrating that: (1) small changes in the MSSM input parameters can lead to 
significant topological changes in the pattern observed on the wedgebox plot; (2) the signal 
and background events often have markedly different distribution patterns on the wedgebox 
plot, pointing toward the possibility of further purifying cuts (perhaps in conjunction with 
extra information garnered from other studies or additional assumptions to clarify of what 
one is obtaining a purer sample) such as the example presented for plot (a) of Fig. 10; and 
(3) the composition of the A° — > Ai + E^'^ signal, that is, what percentages are due to 
H^, X^x'j for different i and j, may be ascertained to some level. The basic topological 

features of the wedgebox plot provide strong, often easily interpreted, leads as to which modes 
are the dominant contributors. The locations of the edges of such features on the wedgebox 
plot also provide information about the sparticle spectrum. The densities of event points in 
each component of wedgebox checkerboard can also be used to distinguish wedgebox plots 
with the same topological features/edges, such as, for instance, telling a wedgebox plot with 
a 2-3 wedge and a 2-2 box from one with only a 2-3 wedge. Further, these point density 
distributions may be used to reconstruct information about the relative production rates of 
the different A° XiX'j processes, though extracting such 'dynamical' information may 
well be far more complicated than is the task of extracting 'kinematical' information about 
the sparticle spectrum from the locations of the edges. All of this is further complicated by 
the remaining background events, and a more holistic study looking at both the Higgs boson 
produced signal and the MSSM backgrounds together may be most appropriate [48]. 
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Note 



Motivated in part by the earlier archival submission of this work, a similar analysis was 
eventually carried out by a member of ATLAS [54], also aiming at mapping out MSSM Higgs 
boson discovery regions via H^, — > XiX'j decays. Results of this ATLAS analysis are 
essentially consistent with those presented herein, though the actual shapes of the discovery 
regions obtained differ somewhat. These differences are in part attributable to adopting 
different selection criteria and employing different simulation tools. Of particular note are the 
ti and hbZ^^*^ backgrounds which are quite significant in the case of the ATLAS analysis but 
yield no background events in this study^^. This is mainly due to the more stringent lepton 
isolation criteria adopted for this study which are very effective at removing leptons produced 
in these two would-be background processes from S-mesons decays. The restrictions on E^, 
which are absent from [54], also aid in removing residual background events. 
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